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FOREWORD 


The  theory  of  optimal  search,  as  usually  presented,  treats  the  development  of  a  search  effort 
allo<  ar  ion  over  a  space  X  which  is  optimal  in  some  sense,  given  an  underlying  target  position 
distribution.  In  other  words,  when  the  location  of  a  target  is  not  precisely  known  and  when  the 
detection  characteristics  of  the  search  mechanism  are  well-defined,  one  can  prescribe  a  search  effort 
function  of  both  space  and  time  to  achieve  a  practical  MOE  (measure  of  effectiveness).  The  best 
measure  of  effectiveness  is  usually  considered  to  be  the  probability  of  detection,  or  target  acquisition, 
under  a  time  budget  constraint. 

Optimal  search  theory,  surprisingly,  usually  does  not  embrace  in  a  rigorous,  methodical  fashion 
the  generation  of  vehicle  path  processes  which  reproduce,  to  some  level  of  approximation,  the 
optimal  allocation  of  effort.  Rather,  more  often  than  not,  heuristic  search  mechanisms,  such  as 
parallel  sweeping  or  radial  tracking,  are  designed  to  replicate  the  effort  distribution.  The  purpose 
of  the  report  is  primarily  to  examine,  from  first  principles,  the  underlying  variational  problems 
and  to  obtain,  if  possible,  the  best  solutions  satisfying  reasonable  cost  criteria.  We  show  that 
it  is  possible  to  generate  sensible  paths  when  one  knows  that  the  target  is  stationary.  Also,  those 
conditions  under  which  the  best  solution  can  be  approximated  by  a  logarithmic  spiral  are  presented. 
Finally,  in  the  case  of  a  moving  target,  we  present  a  stochastic  vehicle  search  pattern ,  whose  expected 
path,  for  example,  could  be  followed  by  a  search  vehicle. 

The  present  document  treats  the  terminal  search  phase  of  an  overall  engagement  strategy 
synthesis  model  for  acquistion  of  an  enemy  target.  The  work  was  sponsored  by  the  Naval  Surface 
Warfare  Center  (NSWC)  under  an  independent  research  grant,  the  principal  investigator  for  the 
overall  effort  being  Dr.  B.  C.  Meyers  (N35). 
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CHAPTER  1 
INTRODUCTION 

In  a  standard  work  on  search  and  detection,1  the  following  comments  are  found:  “It  is  natural 
to  expect  that  a  subject  called  ‘Search  Theory’  will  be  mainly  preoccupied  with  assessing  the 
consequences  of  following  a  particular  search  path,  and  perhaps  even  with  discovering  paths  that 
are  in  some  sense  optimal  for  the  searcher.  This  is  not  the  case-such  problems  are  rare  in  Search 
Theory.. .There  are  two  good  reasons  for  this:  One  is  that  most  searchers  would  not  benefit  from 
knowing  the  optimal  path  unless  it  happened  to  be  easy  to  follow,  and  the  other  is  that  a  ‘path’  is 
not  a  convenient  mathematical  object  in  most  cases.”  This  kind  of  philosophy  may  have  some  merit, 
especially  with  regard  to  the  statement  that  most  searchers  would  not  benefit  from  knowledge  of 
an  optimal,  complicated  path,  which  would,  in  all  likelihood,  not  be  conveniently  traversable  by  a 
searcher.  On  the  other  hand,  it  can  be  argued  that  an  understanding  of  the  optimal  track  is  at 
least  useful  insofar  as  it  provides  a  benchmark  against  which  the  performance  of  a  “simple  track 
process”  can  be  evaluated.  Furthermore,  one  might  be  able  to  approximate  the  complicated  path 
by  a  reasonably  simple  one,  say,  by  a  chordal  or  polygonal  line  approximation  to  the  optimal  track. 
As  to  the  comment  that  an  optimal  path  is  usually  not  a  convenient  mathematical  object,  hopefully 
we  shall  be  able  in  this  document  to  provide  convincing  evidence  to  the  contrary. 

Indeed,  in  this  report,  we  shall  utilize  optimal  control  theory  to  develop  optimal  path  solutions 
for  the  problem  of  acquistion  of  a  stationary  target.  It  will  be  seen  that  such  theory  can  be  invoked 
to  yield  a  “convenient  mathematical  object,”  albeit  someone  might  argue  that  he  would  not  consider 
the  path  to  be  practical.  The  development  of  this  part  of  the  work  is  due  to  the  second  author 
of  this  report,  J.  W.  Wingate,  whose  expertise  is  in  the  field  of  optimal  control.  The  first  author, 
E.  A.  Cohen,  Jr.,  analyzed  the  logarithmic  spiral  search  and  compared  that  path  with  the  optimal 
one  insofar  as  asymptotic  behavior  is  concerned.  The  dynamic  model  was  also  developed  by  the 
first  author. 

It  is  our  hope  that  the  methodology  presented  here  will  serve  two  purposes:  (1)  fill  what  we 
consider  to  be  an  obvious  hiatus  in  the  search  theory  literature  and  (2)  provide  a  theoretical  founda¬ 
tion  for  better  practiced,  computational  studies.  Heretofore,  authors  2'3'4  seem  to  have  emphasized 
various  practical  search  mechanisms;  they  apparently  have  not  addressed  the  issue  as  to  behavior 
of  their  search  mechanisms  in  comparison  with  some  standard. 
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CHAPTER  2 

OPTIMAL  CONTROL  SOLUTIONS  FOR  A  STATIONARY  GAUSSIAN  TARGET 

GENERAL  THEORY 

Before  addressing  the  particular  problems  of  concern  to  us,  it  is  useful  to  review  some  of  the 
theoretical  optimal  control  theory  concepts  pertinent  to  search.  To  that  end  let  us  consider  the 
following  special  case  of  the  general  control  problem  of  Lagrange:  5 

Minimize  /  L(t,x(t),u(t))dt  subject  to  the  constraints  on 
Jo 

x(<),u(<)  that 

x(t)  =  f(f,x(t),u(0),0  <t<T,  (2-1) 

u (t)  €  U ,  (2-2) 

x(0)  =  a;  g\x(T))  =  0,\<i<m,  (2-3) 

where  U  is  a  set  of  vector  functions  u(t)  (called  controls)  and  the  g',  1  <  i  <  m,  are  scalar  functional 
forms  to  be  equated  to  zero  at  the  terminal  time  T  (terminal  constraints).  In  this  special  case 
of  the  control  problem  of  Lagrange,  the  initial  vector  a  is  given  and  the  interval  [0,T]  is  fixed. 
Furthermore,  the  vector  x  belongs  to  an  n-dimensional  Euclidean  space  Rn,  and  u  belongs  to  a 
k-dimensional  Euclidean  space  Rk.  The  functions  L  and  f  are  to  be  continuous  in  t,x,  and  u  and  to 
have  continuous  first  partial  derivatives  L,,LXl  f,,  and  fr.  That  is,  L  and  f  are  required  to  be  smooth 
functions.  In  addition,  the  functions  g *  are  to  be  smooth.  In  contrast,  the  control  u(<)  will  generally 
be  assumed  to  be  piecewise  continuous. 

For  this  special  case,  the  Maximum  Principle  of  Optimal  Control  Theory  takes  the  following 

form: 

Define  the  function  H  of  (t,x,u,  p)  by 
H(t,  x,  u,  p)  =  p  •  f(t,  x,  u)  -  A0L(<,  x,  u), 

where  A0  is  a  scalar  constant.  The  function  H  is  known  as  a  Hamiltonian.  If 

Xo  :  xo(t),uo(<),0  <  t  <  T, 

is  a  solution  to  our  minimization  problem  stated  at  the  outset,  there  exist  multipliers  A;,0  <  j  <  m. 
not  all  zero,  and  an  absolutely  continuous  function  p(t)  (understood  componentwise)  such  that 

A0  >  0;  (2-4) 
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the  Euler-Lagrange  equations 

xo(<)  =  Hp(t,xo(t),u0(t),p(t)) 

p(t)  =  -Hx(t,  x0 (<),  u0(t),  p(t)),  (2-5) 

0  <  t  <T,  hold  (the  dot  denoting  differentiation);  the  Weierstrass  condition: 


tf(t,xo(t),uo(0,p(t))  =  max  H(t,xo(t),\i,p(t)),0  <t<T  (2-6) 

obtains;  the  equation 

^■(i/(t,xo(t),uo(t),p(<))  =  Ht(t,  xo(t),  u0(t),  p(t)),  0  <t  <T  (2-7) 

holds  (d/dt  denoting  a  total  derivative  and  Ht  representing  a  partial  derivative  with  respect  to 
<-x,  u,  and  p  being  held  fixed);  xo  satisfies  the  imposed  boundary  conditions  given  by  Equation  (2- 
3);  and  p  satisfies  the  natural  boundary  (or  transversality  conditions: 

m 

EA^xW7'))  +  p(T)  =  0.  (2-8) 

i=i 

If  A0  >  0),  it  can  be  normalized  to  a  convenient,  positive  constant  (say  1).  When  A0  =  0,  the  solution 
is  said  to  be  abnormal.  Such  abnormality  usually  arises  when  too  many  constraints  are  imposed. 

In  practice,  Differential  Equations  (2-5)  will  often  be  complex  and  would  need  to  be  solved 
numerically  if  a  direct  approach  were  necessary.  If  the  solution  to  this  system  can  be  inferred  by 
other  means,  say  through  qualitative  properties  of  the  system  (sign  properties  of  right-hand  side, 
for  example)  and  through  characteristics  of  the  Weierstrass  condition,  namely,  Equation  (2-6),  then 
we  are  able  to  ascertain  the  answer  indirectly.  We  shall  see  the  advantages  of  such  artifices  in  the 
sequel. 


THE  OPTIMAL  SEARCHER  MOTION  PROBLEM 

Let  us  now  investigate  the  details  of  the  optimal  searcher  motion  problem  from  the  perspective 
of  optimal  control  theory.  We  shall  develop  optimal  solutions  for  two  practical  situations:  (A)  The 
target  is  stationary  in  two-space  (R2),  and  its  location  is  governed  by  a  circular  Gaussian  Law: 


p(x)  = 


(2-9) 


and  (B)  the  target  is  stationary  in  two-space,  and  its  position  follows  an  elliptic  Gaussian  Law: 

p(x,=  2T^'*p(~K^+f))-  (2'10) 

In  Case  (B),  we  can  always  orient  the  axis  of  our  coordinate  system  (in  principal  directions)  so  as 
to  eliminate  correlation  effects.  We  shall  develop  our  solution  for  Case  (A)  first,  utilizing  the  theory 
of  Lagrange  just  described.  For  (B),  a  different  approach,  subsuming  (A),  will  be  presented. 
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For  Case  (A),  let  r(<)  =  (xi(t),  x2(t))  represent  the  searcher’s  position  as  a  function  of  time,  and 
let  ur  suppose  that  the  searcher  is  to  follow  a  constant  speed  trajectory  which  moves  outward  from 
the  origin  (the  mode  of  the  distribution).  Then  the  searcher  motion  is  described  by 

r(0  =  v(<).(r  G  R2), 

where  ||v(T)||  =  v,  and  r  v  >  0.  The  latter  condition  follows  from  the  requirement  that  ||r||  is 
nondecreasing,  since  d(\\r\\2)/dt  =  c/(r  r)/dt  =  2r  •  v  >  0. 

Two  issues  pertinent  to  Problem  (A)  need  to  be  addressed:  (1)  an  optimal  allocation  of 
effort  corresponding  to  the  target  distribution  and  (2)  a  cost  effectiveness  criterion  to  measure  the 
efficiency  of  any  given  search  path  in  matching  the  optimal  allocation.  Point  ( 1 )  has  been  resolved 
in  the  literature.6  We  shall  first  review  that  argument. 

The  target  distribution  function  for  Problem  (A)  has  been  given  by  Equation  (2-9).  One  also 
needs  a  detection  law,  which  gives  the  probability  of  acquisition  of  a  target  known  to  be  located  at 
some  position  x  as  a  function  of  x  and  an  effort  density  z.  The  simplest  law  of  exponential  type  is 
afforded  by 

b(x,z)  =  1  -  e~az,  x  €  X,a  >  0,  z  >  0,  (2-11) 

where  a  is  some  known  quantity.  Note  that  b(x,z),  as  given  by  Equation  (2-11),  does  not  depend 
explicitly  upon  x.  However,  there  may  well  be  an  implicit  dependence  in  the  sense  that  z  could  be 
a  function  of  x.  In  fact,  one  would  generally  expect  that  the  latter  functional  coupling  would  be 
present.  Equation  (2-11)  is  certainly  reminiscent  of  the  Koopman  law  of  random  search,2  but  it 
differs  conceptually  in  the  sense  that  Koopman ’s  Law  gives  an  unconditional  probability  of  detection. 
whereas  Equation  (2-11)  is  definitely  a  conditional  one.  It  is  important  to  realize  that  the  optimal 
allocation  which  we  derive  now  is  predicated  on  the  kind  of  detection  law  assumed.  If  Equation 
(2-11)  turned  out  to  be  an  inadequate  description,  one  would  have  to  substitute  for  it  another 
detection  law  and  consequently  derive  another  allocation.  Conceptually,  the  methodology  which 
will  be  presented  should  not  generally  be  affected;  only  the  optimal  effort  distribution  need  be 
modified. 

The  quantity  z  of  Equation  (2-11)  is  to  be  an  effort  density,  i.e. ,  in  our  case,  2  is  to  be  given 
by  2  =  /(x)  such  that 

I  f(x)dx  —  A,,  (2-12) 

Jx 

where  A,  is  the  total  area  to  be  searched,  using  a  sensor  of  width  H'.  Explicitly,  A,  =  It'/,,  where  L 
is  the  length  of  the  search  path.  Thus  z  is  dimensionless  and  is  indeed  an  areal  density  of  search 
in  the  truest  sense  of  the  term.  Now,  if  one  uses  an  effort  density  /(x )  and  any  giv^n  conditional 
detection  law  6(x,/(x))  for  that  effort  density,  one  can  express  the  unconditional  chance  of  detection 

by 

P(f)=  f  p(x)b(xj{x))dx. 

Jx 


(2-13) 
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so  that  P(f)  is  obviously  a  functional  (a  mapping  of  functions  to  the  real  numbers).  We  want  to 
obtain  a  function  /*(x)  which  maximizes  P(f)  under  the  constraint  of  Equation  (2-12).  That  is,  we 
want  an  allocation  of  resources  which  maximizes  the  acquisition  probability.  Note  that  a  “cookie- 
cutter  effect”,  in  which  6(x,z)  ~  1  for  a  negligible  effort  expenditure  z,  is  achieved  by  choosing  a  to 
be  a  large  positive  number. 

The  rest  of  the  argument  for  Issue  (1)  above  is  furnished,  following  Stone,6  as  follows:  for  a 
fixed  A  >  0,  one  maximizes  the  pointwise  Lagrangian  with  respect  to  z,  given  x: 


L(x,z)  =p(x)(  1  -e  ai)  -  Az. 


One  finds  that  the  optimal  solution  Z*  =  fl  is 


**  =  /a(x)  =  -0nMx)/A]}+, 

a 


where 


M8' 
1 0, 


Theory  presented  by  Stone  then  indicates  that  fl  is  optimal  over  all  functions  /  satisfying 


J^f(x)dx  <  J*fl(x)dx. 


(2-14) 


Therefore,  if  one  finds  a  A  >  0  such  that  Equation  (2-12)  holds,  he  will  have  the  optimal  allo:ation. 
The  remainder  of  the  argument  follows  that  of  Stone,6  there  worked  out  for  the  case  a  =  1.  We 
have,  first  of  all,  the  following  optimal  allocation,  given  in  polar  coordinates  (r,0)  : 


f-  —  ln(2x<r2Ai)  —  ,r2  <  -2<r2  ln^ir^A!), 

fxt(r,6)  =  «[  2<r2J 

1  0  ,r2  >  — 2<r2  ln(2jra2Ai), 

where  A]  =  A/a.  Integrating  the  cost  function  over  the  space  X,  as  Stone  does,  we  find  that 

C(/a,)  =  xo-2[ln(2x(T2A1)]2/a. 

Equating  A,  to  the  right  hand  side  of  Equation  (2-16)  and  solving  for  A,,  we  have 


(2-15) 


(2-16) 


•'‘  =  5^'xp  f"  ($) 


Substituting  Equation  (2-17)  for  Ai  into  Equation  (2-15),  we  find  that 


);[(£)*-£]  (3)' 

1° 


(2-18) 
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Note  that  the  argument  presented  so  far  does  not  require  a  constant  speed  search.  If  we  do  assume 
that  the  searcher  travels  at  constant  speed  v,  then  A,  =  WvT,  where  T  is  the  total  time  devoted  to 
the  search.  Observe  that,  at  any  point  (r,0),/J(r,0)  =  0(  1/ai)  as  a  tends  to  -t-oo,  i.e.,  fl(r,e)  tends  to 
zero  as  a  tends  to  +oo;  it  does  so  like  a~K  Thus,  as  the  detection  law  approaches  that  of  a  cookie 
cutter  (for  any  effort  z ,  the  detection  probability  is  1,  provided  the  target  is  at  x),  the  effort  needed 
at  any  point  tends  to  zero,  but,  at  the  same  time,  the  effort  is  supported  over  a  disk  whose  radius 
tends  to  infinity  like  aE  Clearly,  the  searcher,  under  these  circumstances,  should  “move  out  from 
the  origin  more  and  more  rapidly,”  i.e.,  he  should  spend  very  little  time  in  any  given  area  trying  to 
locate  a  target.  The  mathematics  seems  to  reflect  this  intuitive  line  of  reasoning. 

Suppose  that,  instead  of  assignation  of  effort  over  the  plane,  one  decides  a  priori  to  restrict 
his  attention  to  some  disk  D  centered  on  the  origin.  Now  P(D ),  the  chance  that  the  target  lies  in 
D,  is  given  by 

p(0)=2 hS,  l  exp(i£)r<,r'" 

=  1  -  “p  (-§,)  <m9> 
Thus,  from  Equation  (2-19),  if  P  =  P(D)  is  given  in  advance,  we  can  easily  find  R.  In  fact,  R  = 
<r[2  ln(  1/(1  -  P))]E  Now,  assuming  that  the  target  belongs  to  D,  the  target  distribution  becomes  a 
conditional  one,  namely, 


Pd(x) 


1  (  x\  +  xl\ 

2  n^P(D)eXp{-^-) 
0 


,x  6  D 
,xfO. 


(2-20) 


With  this  new  target  distribution,  one  now  attempts  to  find  the  optimal  allocation  of  effort.  We 
have 


zd  =  /a,d(*)  =  ^{ln[apn(x)/A]}+,x  €  D, 


so  that,  if  R2(P,  Ai ,  a)  =  min(E2,-2<r2  ln(27r<72A1P)), 
/;,nM)  =  (  « 


.2  ' I 


-ln(2x<r2A XP(D)  -  ^ 


i  r2  <  R*(P,\u<r) 
l  0  ,  r2  >  R2(P,  A, ,  a). 

Thus,  if  we  define  A2  =  XiP(D),  we  find  that  the  cost  over  the  disk  is 


c(/;,dm))=  l 


1—  ln(2ir<r2A2)  - 


4  <T2 


] 


<  Aj  < 


1  -  P 
2 Kcr2  ' 


v  ^~[ln(2x<72A2)]2 


1  -  P 


2  X<72  —  ^2  —  2t Tff2 

The  first  line  of  Equation  (2-23)  corresponds  to  the  case  where  R?(P,X\,<r)  =  R}.  In  that 
event  D  is  nonzero  over  the  entire  disk  p  <  R.  In  fact,  0  <  A2  <  (1  -  P)/ 2-na2  if  and  only  if 


1 


(2-21) 


(2-22) 


(2-23) 
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R2  <  — 2<r2  ln(27rcr2A2).  Equating  the  right  side  of  line  1  of  Equation  (2-23)  to  A,  =  WvT ,  one  finds 
that  A2  =  exp  Inserting  this  relation  for  A2  into  the  inequality  0  <  A2  <  (l-P)/2*a7, 

one  finds  that  A,  >  it a2  ln2(l  -  P)/a.  Also,  the  allocation  density  can  be  found  in  terms  of  WvT  by 
using  in  line  1  of  Equation  (2-23)  the  same  equality  for  A2.  It  follows  that 

fx,Dir^)  =  ~  tJt~  )'r  <  R,WvT>  ir<72(ln2(l  -  P))/a.  (2-24) 

In  addition,  we  would  require  that  R  <  vT ,  so  that  a  searcher  could  reach  the  circle  p  =  R  at  time  T. 

On  the  other  hand,  consider  the  second  line  of  Equation  (2-23).  Then  we  find,  of  course, 
that  0  <  A,  <  ttct2  ln2(l  -  P)/a,  corresponding  to  <  A2  <  Now  A2  =  exp  (-g-( 

Inserting  this  equality  into  the  allocation  density,  we  have,  when  WvT  <  ?r<r2(ln2(l  -  P))/a  and 
R\  S  —  2<t2  ln(2ircr2A2)  =  2<r(h^)r, 


fx,oW) 


(  1  fWvTa\ *  r2 

w  )  2  <r2J 


,  r2  <  R\  <  R2 
,R\  <  r2  <  R2. 


(2-25) 


Equation  (2-25)  has  the  same  form  as  Equation  (2-18),  the  only  difference  being  that,  instead 
of  having  a  location  distribution  over  the  entire  plane,  we  have  one  restricted  to  the  disk  r  <  R. 
When  the  effort  to  be  expended  does  not  exceed  the  threshold  ir<r2(ln2(l  -  P))/a,  we  have  a  situation 
similar  to  that  encountered  when  the  size  of  the  disk  is  unrestricted.  Any  theory  developed  for  that 
case  would  then  be  applicable  for  Equation  (2-25)  also.  The  situation  for  Equation  (2-24)  is  quite 
interesting;  however,  further  details  concerning  it  will  be  deferred  for  the  moment. 

Let  us  now  return  to  Problem  (A).  The  polar  form  presented  in  Equation  (2-18)  will  be  quite 
useful  to  us  in  the  development  of  a  cost  effectiveness  model,  which  is  the  goal  of  Step  (2).  Of  course, 
the  radially  symmetric  nature  of  Problem  (A)  introduces  an  element  of  simplicity  upon  which  we 
can  certainly  capitalize.  Note,  for  example,  that  Equation  (2-18)  includes  no  9  dependence.  Even  so, 
generation  of  a  cost  criterion  which  is  ideally  suited  to  the  task  of  path  effectiveness  assessment  is  a 
complex  issue.  In  some  sense,  one  would  like  to  match  the  allocation  density  prescibed  by  Equation 
(2-18)  with  the  allocation  density  which  evolves  from  any  particular  choice  of  path.  However,  as 
Figure  2-1  indicates,  the  latter  density  is  a  very  complicated  entity  because  of  the  complex  fashion 
in  which  any  particular  spatial  point  is  generally  covered  by  search  area,  even  for  a  path  with 
increasing  radius  vector  in  time.  Note,  in  Figure  2-1,  that  calculation  of  the  actual  search  density 
would  involve  contributions  from  several  strips  of  width  W  and  that,  the  closer  the  spiral  arms  are 
to  one  another,  the  harder  such  a  calculation  would  become. 

Clearly  we  need  a  cost  criterion  based  on  a  simpler,  but  still  meaningful,  concept.  Therefore, 
instead  of  concerning  ourselves  with  the  allocation  density,  let  us  attempt  to  match  a  measure 
of  cumulative  allocation.  Since  Problem  (A)  exhibits  radial  symmetry,  one  is  inclined  to  match 
cumulative  allocations  within  disks  centered  on  the  origin.  That  is  to  say,  suppose  that  a  curve 
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FIGURE  2-1.  SEARCH  PATH  DENSITY  ASSOCIATED  WITH 
TYPICAL  POINT  x 

such  as  that  displayed  in  Figure  2-1  is  expressed  in  polar  form  as  a  function  of  time,  given  by 
(p(t),0(t)),  and  that  we  agree  to  compute  the  total  area  searched  using  such  a  path  up  to  time  t  -  t0. 
For  a  constant  speed  of  search  v,  that  area  would  just  be  Wvt0.  The  cumulative  effort  function 
to  be  matched  within  the  disk  of  radius  p(t0)  is  obtained  by  integration  of  the  right  hand  side  of 
Equation  (2-18)  over  that  disk.  We  denote  the  latter  quantity  by  E(p{t0)),  where  p(io)  =  ||r(<0)||-  To 
simplify  matters,  we  agree  to  ignore  search  area  which  intersects  the  disk  {r  :  ||r||  <  p(t0)}  for  time 
t  >  to-  Hopefully,  that  contribution  is  relatively  small.  Thus,  we  are  led  to  the  error  (cost)  criterion 
of  mean-square  type: 

C(r(<))  =  f  (c(t,p(t)))2  dt,  c(t,p)  =  Wvi  -  E(p),  (2-26) 

Jo 

which  is  at  least  a  tractable  quantity.  Criterion  (2-26)  does  suffer  from  one  difficulty,  viz.,  that 
peculiar  oscillatory  paths  r (t)  could  arise  which  optimize  it  despite  the  fact  that  such  tracks  are 
unrealistic  candidates.  To  alleviate  that  situation,  one  might  require  that  both  p{i)  and  9(t)  be 
strictly  increasing  functions  or  that  p(t)  be  strictly  increasing  and  0(t)  be  strictly  decreasing. 
Integration  of  Equation  (2-18)  over  a  disk  of  radius  p(t)  leads  to  the  quantity 


BMO)  =  (#(T)  -  ? 22). 


(2-27) 


where 


R2(T) 
2  a2 


faWvTY 

V  ™’2  / 
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We  demand  that  our  path  start  from  the  origin  and  intersect  the  circle  ||r||  =  R(T)  at  time  t  =  T. 
That  is,  we  have  the  boundary  conditions  r(0)  =0;||r(T)||  =  R(T).  The  quantities  W,v,T,a,  and  R  are, 
for  the  present  purpose,  to  be  treated  as  given  positive  numbers. 

The  circular  symmetry  of  problem  (A)  allows  us  to  reduce  a  2-D  problem  to  a  1-D  one  in 
terms  of  the  radial  variable  p  =  ||r||.  Anticipating  future  developments,  after  production  of  the  radial 
solution  part,  the  azimuthal  motion  can  be  determined  by  choosing  the  azimuthal  speed  such  that 
the  total  speed  is  v.  Note  that  v 2  =  i2(t)+y2(t)  =  p2(t)+p2(t)92(t)  (where  a  dot  denotes  differentiation), 
so  that  |/?(t)|  <  t>,  a  natural  restriction  on  the  radial  rate.  Therefore,  the  radial  motion  satisfies 

p(i)  =  u, 

where  0  <  u  <  v.  Also,  p(0)  =  0  and  p(T)  =  R(T).  The  cost  functional  to  be  considered  is  that  of 
Equation  (2-26).  Let  us  set  c(t,p(t))  =  0  for  every  t.  Appealing  to  Equation  (2-27),  one  then  finds 
that 

p(t)  =  fi(T)(l-(l-t/T)i)i.  (2-28) 

Shown  in  Figure  2-2  is  the  curve  p  -  p(t),  together  with  those  regions  for  which  c  <  0  and  c  >  0. 
Since  p(t)  =  +oo  at  both  t  =  0  and  t  -  T,  it  is  clear  that  our  speed  constraint  on  the  radial  motion 
will  not  allow  us  to  secure  a  path  for  which  c(t,p(t))  is  always  zero,  i.e.,  p(t)  cannot  be  followed  near 
the  endpoints. 


T 


FIGURE  2-2.  CURVE  p  =  p{t)  AND  REGIONS  c  <  0  AND  c  >  0 

We  now  develop  the  necessary  conditions  for  optimality,  using  the  optimal  control  theory 
results  of  Inequality  (2-4)  and  Equations  (2-5)  to  (2-8).  The  Hamiltonian  in  our  case  is 

H  =  p(t)u  -  A0[c(t,*>)]2. 
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We  are  looking  for  an  optimal  pair  po(t),  «o(0-  Hamiltonian  theory  shows  that  one  can  find  an 
absolutely  continous  function  p{*)(0  <t  <  T)  and  a  nonnegative  multiplier  A0  such  that 

(Ao,p(0)#(0,0)  (0  <t<<T) 

and  such  that  the  Euler- Lagrange  equations  hold: 


Po(i)  =  «o(0 

p(t)  =  — 2A0c(t,  Po(t))E'(po(t)), 


0  <  t  <  T. 


(2-29a) 

(2-29b) 


(Both  the  dot  and  the  prime  denote  differentiation  operations.  Throughout  the  text,  we  shall  utilize 
both  notations  interchangeably.)  In  addition,  we  have  the  boundary  conditions 


Po(0)  -  0;po(T)  =  R 

mentioned  previously  and  the  Weierstrass  condition  of  Equation  (2-6): 


(2-30) 


p(t)u0(t)  -  Ao[c(t,p0(<))!2  =  I™**  (p(<)“-  A0[c(bpo(<))]2). 


Thus 


«o  (<) 


=  {V' 
l  0, 


(2-31) 


P(t )  >  0 

p(t)  <  0, 0  <  <  <  T, 

no  information  being  provided  about  u0(<)  when  p(t)  =  0.  The  condition  where  p(t )  =  0  is  known  as 
a  singular  control  situation. 

The  multiplier  A0  must  be  nonnegative.  Suppose  that  A0  =  0.  Then,  of  course,  p(t)  can  never 
vanish,  and  its  continuity  implies  that  it  has  the  same  sign  for  all  t  in  [0,T].  If  p(t)  were  negative  for 
all  t,  Equation  (2-31)  would  indicate  that  u0(t)  =  0  everywhere.  Then  Equation  (2-29a)  leads  to  the 
conclusion  that  p0(t)  is  constant  for  all  t.  Since  po(0)  =  0,  indeed  po(t)  =  0  for  all  t.  The  trajectory 
would  be  reduced  to  a  single  point,  and  obviously  a  one-point  trajectory  is  not  a  solution  to  our 
problem.  Thus  the  sign  of  p(t)  is  positive  for  all  t.  Now  Equation  (2-31)  implies  that  uo(<)  =  v  for  all 
t  €  (0,T).  Appealing  again  to  Equation  (2-29a)  and  Equation  (2-29b),  together  with  the  fact  that 
po(0)  =  0,  one  has  purely  radial  motion  p  =  vt, 6  =constant.  In  order  that  we  meet  the  boundary 
condition  p(T)  =  R{T),  we  must  have  R(T)  =  vT.  Now  this  is  a  tenable  motion  (and  indeed  the  only 
one  possible)  when,  in  fact,  R(T)  =  vT.  One  must  admit  that  this  instance  is  trivial;  so,  without 
loss  of  generality,  let  us  assume  that  R  <  vT. 

Assume  then  that  A0  >  0.  As  mentioned  previously,  we  may  normalize  it  to  be  unity.  Let 
us  suppose  that  p(T)  <  0.  Then,  since  p(t)  is  continuous,  p(t)  <  0  within  some  c  neighborhood  of 
t  =  T,  say  on  (T-  e,  T).  Condition  (2-31)  shows  that  u0(<)  =  0  on  the  c-interval,  and  Equation  (2-29a) 
says  that  p0(t)  is  constant  on  the  same  interval.  The  boundary  condition  p0(T)  =  R  implies  that 
po(0  =  R,t  €  (T-  t,T).  We  assert  that  p0(t)  =  R,0  <  t  <  T.  Suppose,  to  the  contrary,  that  there  exists 
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a  t,  say  t  =  <1,  where  ti  <T  -  e,  and  that  /jo(<i)  <  R-  By  continuity  of  po(t),  we  may  assume,  without 
loss  of  generality,  that  (*,p0(<))  lies  in  the  region  c  <  0  of  Figure  2-2  for  t  >  *i.  Appealing  to  Equation 
(2-29b)  and  the  fact  that  E'(p0(t))  >  0,<  >  <i,  one  finds  that  p(t)  >  0 ,<  >  <i.  Therefore,  p(t)  <  0  in 
Condition  (2-31)  shows  that  u0(t)  =  0  in  [tuT\.  Thus,  p0(t)  is  constant  in  [fi.T’J.  Since  po(t)  =  R  when 
t  €  (T  -  e,T),p0(t)  =  R  in  [tltT].  We  have  a  contradiction.  It  follows  that  po(t)  =  R,0  <  t  <  T.  From 
Equation  (2-29a)  and  Equation  (2-29b),  since  E'(p0(t))  =  E'(R )  =  0,p(t)  =  0  for  all  t,  and  so  p(t)  =  p(T). 
We  have  po(t)  =  R,  u0(t)  =  0,  p(t)  =  p(T)  for  0  <  t  <  T.  Such  a  state  of  affairs  is  in  conflict  with  the  fact 
that  po(0)  =  0.  Suppose  that  p(T)  =  0.  We  claim  that  again  there  exists  an  e  >  0  such  that  u0(<)  =  0 
when  t  belongs  to  ( T  -  c,  T).  Suppose,  to  the  contrary,  that,  for  every  c  >  0,  there  exists  t<  such 
that  u0(t()  >  0.  Now  we  know  that  there  exists  a  6  >  0  such  that  c(t,p0(t))  <  0  in  (T  -  6,T).  We  see 
that  this  is  true  by  checking  Figure  2-2  and  noting  that,  unless  the  curve  p  =  po(t)  lies  in  the  region 
c  <  0  for  t  sufficiently  close  to  T,  the  radial  rate  requirement  would  be  violated.  In  fact,  the  only 
way  p  =  p0(t)  could  lie  completely  in  c  >  0  for  t  close  to  T  is  that  a>(T)  =  +oo.  Such  a  condition  is  a 
clear  violation  of  physical  requirements.  Also,  if  c  -  0  at  a  point  t0  sufficiently  near  t  =  T,  then  since 
p(T)  =  +oo,  the  slope  of  the  line  joining  (t0,p(<0))  to  ( T,p(T ))  must  exceed  v.  The  law  of  the  mean 
for  derivatives  then  implies  the  existence  of  a  point  ti  >  t0  such  that  po(*i)  >  v  (a  clear  violation). 
Suppose  then  that  we  have  t  =  t(  such  that  T  -  6  <  t(  <  T  and  u0(<t)  >  0.  Clearly,  p0{t()  <  R,  for 
otherwise  there  would  exist  a  t  >  tc  for  which  po(t)  >  R  (by  the  definition  of  a  derivative).  Thus, 
E'(po(tc))  >  0,c(tc,  p0(t())  <  0,  and  A0  >  0,  so  that  p(te)  >  0.  Now,  appealing  to  Equation  (2-29a)  and 
Equation  (2-29b)  again,  p(t)  >  0  in  [tf,T].  Therefore,  p(i)  is  monotone  increasing  in  that  interval, 
and,  since  p(T)  =  0,p(<)  <  0  in  [<flTj.  Suppose  that  p(te)  =  0.  Since  p(te)  >  0,  there  would  exist  a  t  >  t( 
such  that  p(t)  >  0,  which  is  impossible.  Hence,  p(tf)  <  0.  From  Equation  (2-31),  u0(t()  =  0,  and  we 
have  a  contradiction  to  our  original  statement.  Therefore,  there  exists  some  f  >  0  such  that,  for 
t  G  [T  -  e,  7],  u0(t)  =  0.  The  rest  of  the  argument  proceeds  as  for  the  case  where  p(T)  <  0.  We  are 
again  led  to  the  fact  that  our  initial  condition  p( 0)  =  0  cannot  be  satisfied.  Finally,  then,  p(T)  must 
be  positive,  and  the  continuity  of  p(t),  together  with  Equation  (2-31),  means  that  u0(<)  =  v  on  the 
final  part  of  the  motion.  As  we  have  shown  previously,  c  <  0  on  the  final  part  of  the  motion,  so  that 
Equation  (2-29a)  and  Equation  (2-29b)  can  be  invoked  once  more  to  show  that  p(i)  >  0  there. 

Remember  that,  as  long  as  p(t)  remains  positive,  u0(<)  =  v,  so  that  motion  is  radial,  and 
thus,  in  Figure  2-2,  po(t)  is  a  straight  line  running  backwards  from  the  point  (T,R(T)).  Thus,  we 
are  within  the  region  c  <  0  until  the  line  po  =  v(t  -  T)  +  R{T)  intersects  p  =  p(t).  We  said  that,  on 
this  part  of  the  motion,  p(t)  >  0.  It  follows  that  p(<)  decreases  as  we  work  backward  along  the  line 
p0  =  vt  k{k  =  R(T)  -  vT)  from  the  point  ( T ,  R(T)).  Radial  motion  will  cease  only  if  p(t)  =  0  for  some 
t.  Let  us  examine  the  first  t  for  which  this  phenomenon  occurs.  (If  p(t)  >  0  for  all  t,  radial  motion 
p0(<)  =  «/<,0  <  t  <  T,  is  the  only  motion,  and  then  our  situation  is  again  the  trivial  scenario  in  which 
R  =  vT.)  We  have  three  possibilities:  (a)  p(7)  =  0  and  c(t,p0(t))  <  0(?  <  t  <  T);  (b)  p(<)  >  0  when 
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the  path  backwards  from  time  T  crosses  c  =  0;  (c)  p(t)  vanishes  just  when  the  backwards  trajectory 
(following  a  straight  line  with  slope  v)  crosses  the  curve  c  =  0. 

First  we  check  Case  (a).  In  this  instance,  from  the  Euler- Lagrange  equations,  p(7)  >  0,  so  that 
p  changes  sign  at  7.  Therefore,  p(t)  <  0  in  some  interval  (7  -  6,7).  From  the  Weierstrass  condition, 
tio(<)  =  0  in  that  same  interval,  and  it  follows  that  p0(t)  =  p0(7)  in  (7-6,7).  It  is  now  possible  to  argue, 
in  a  similar  vein  as  before,  that  u0(t)  =  0  on,  [0,7],  so  that  p0{t)  =  p0(7)  for  t  <  7.  The  initial  condition 
po(0)  =  0  cannot  be  satisfied,  so  that  (a)  is  untenable.  Figure  2-3  is  an  illustration  of  Case  (a). 


FIGURE  2-3.  p(7)  =  0  AND  c(t,p0(t))  <  0,7  <  t  <  T 

For  Case  (b),  we  are  really  in  an  adverse  situation,  for,  from  the  Euler- Lagrange  equations, 
p  <  0  when  c  >  0.  Thus  p(t)  increases  as  we  work  backward  toward  t  =  0  from  the  crossing  point  and 
hence  is  always  positive.  The  Weierstrass  condition  forces  us  to  continue  on  a  radial  course:  and, 
when  R  <  vT,  we  cannot  satisfy  our  initial  condition  po(0)  =  0.  This  situation  is  illustrated  in  Figure 
2-4  and  Figure  2-5.  In  Case  (b),  p(<)  cannot  vanish. 

P 
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P 


FIGURE  2-5.  BEHAVIOR  OF  />(<)  IN  CASE  (b) 

The  remaining  possibility  is  that  p  vanishes  just  when  the  backwards  trajectory  crosses  the 
curve  c  =  0.  The  trajectory  cannot  cross  c  =  0  at  that  point,  since  we  would  be  forced  back  into  a 
situation  like  that  which  occurs  in  Case  (b),  namely,  that  one  where  p  <  0  for  c  >  o.  A  radial  solution 
would  be  the  result,  and  that  possibility  cannot  occur,  since  again  p  would  vanish  at  t  —  T  -  R/v. 
We  must  therefore  follow  the  curve  c  =  0  (so  that  p(t)  =  0  along  it)  until  we  get  to  the  initial  part  of 
the  trajectory,  which  has  the  same  properties  as  that  for  the  backward  portion  from  (T,R{T)).  The 
situation  is  illustrated  in  Figure  2-6  and  Figure  2-7.  The  solution  has  three  component  parts.  From 
0  up  to  some  time  ti,u0(t)  =  v  and  c  >  0;  from  tl  to  t2,  u0{t)  =  KO.c  =  0,  and  p0(t)  =  p(<);  from  t2  to 
T, uo(t)  =  v  and  c  <  0.  In  [0,ti]  p  is  positive;  on  [<i,<2]  it  vanishes  (singular  control  problem);  on  [t2,T] 
it  is  again  positive.  On  the  singular  portion,  «0(<)  <  v,  so  that  azimuthal  motion  is  required  to  bring 
the  speed  up  to  v.  On  the  nonsingular  portions,  motion  is  radial.  The  switching  times  1 ,  and  t2  are 
easily  found,  since  we  know  the  equations  for  p  and  for  the  straight-line  portions  of  p0.  We  have 

vt1  =  R(l-(l-t1/T)i)i  (2-32) 

R-v(T-h)  =  R(  1  -  (1  -  ’  .  (2-33) 

To  solve  Equation  (2-32),  we  first  define  K  =  R/vT, 0  <  K  <  1.  Since  vI<T  =  <  1<T.  In  fact,  from 

Figure  2-8,  one  sees  that  ti  <  KT.  Let  «i  =  ti/KT  =  vti/R.  Then,  from  Equation  (2-32), 

•i  =  (l-Cl-Fsi)*)*,  0  <  s,  <  1, 

so  that  s!  satisfies 

1  -s\  =(!  -/fsi)’,  0  <  si  <  1.  (2-34) 


From  Equation  (2-34),  upon  squaring  both  sides,  we  find  that  si  must  be  a  solution  of  the 
following  equation: 


«(s3  — 2s  +  /O  =  0,0<s<  1. 
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P 


FIGURE  2-6.  BEHAVIOR  OF  p(t)  IN  CASE  (c) 


P 


FIGURE  2-7.  BEHAVIOR  OF  p{t)  IN  CASE  (c) 


Clearly  si  =  0  is  not  the  solution  desired,  for  we  demand  that  0  <  sj  <  1.  Evaluating  the  cubic 
factor  of  the  left  member  of  Equation  (2-35)  at  s  =  -oo,  we  obtain  -oo,  at  s  =  0  we  obtain  K  >  0.  at 
s  =  1  we  have  K  -  1  <  0,  and,  finally,  at  s  =  +oo,  we  get  +oo.  Therefore,  the  cubic  changes  sign  in 
(-oo,0),  (0, 1),  and  (l,oo).  It  follows  that  all  roots  of  the  cubic  are  real  and  simple,  one  zero  occurring 
in  each  of  the  three  aforenamed  intervals.  We  desire  the  root  in  (0,1).  Applying  standard  formulas 
for  solving  cubics,  7  one  has 


=  (2\/6/3)cos  (jfCos"1  (-3*^/8)  +  4w)^  , 


(2-36) 


which  one  can  easily  show  (appealing  to  simple  inequalities)  is  the  smaller  positive  root.  As  is 
customary,  Cos-1  represents  the  principal  arccos  function,  whose  range  is  [0,*].  For  the  solution 
to  Equation  (2-33),  refer  to  Figure  2-9.  It  is  clear  from  the  definition  of  K  that  a  line  with  slope 
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FIGURE  2-8.  FIRST  SWITCHING  TIME  FOR  OPTIMAL 
SEARCHER  MOTION 

v  emanating  from  (T,  R)  must  have  intersected  the  t  axis  at  time  (1  -  I\)T .  We  have  previously 
defined  the  time  t2  to  be  that  time  for  which  the  same  line  intersects  the  curve  p  -  p(0-  Let 
s2  _  (i  -t2/T)/K  =  v(T-t2)/R.  Then,  dividing  Equation  (2-33)  by  R,  one  finds  that 

l-«2=(l-(A*2)*)*,  (2‘37) 

whence 

(Ks2)$  =  2s2  -  s2  (2-38) 

Squaring  both  sides  of  Equation  (2-38)  and  noting  that  s2  ^  0  (equivalently.  t2  ^  T),  we  see  that  s2 
satisfies  the  equation 

s3  —  4s2  +  4  s  —  K  =  0  (2-39) 


P 


FIGURE  2-9.  SECOND  SWITCHING  TIME  FOR  OPTIMAL 
SEARCHER  MOTION 
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The  cubic  polynomial  in  Equation  (2-39)  is  negative  when  s  =  0,  positive  when  s  =  1,  negative 
when  s  =  2,  and  certainly  positive  in  a  neighborhood  of  s  =  -foo.  It  follows  that  all  its  roots  are  real 
and  positive  and  that  there  is  precisely  one  root  in  (0, 1).  Using  a  standard  procedure  for  converting 
a  cubic  to  another  cubic  with  a  missing  quadratic  term,7  let  s  =  u  +  4/3  in  Equation  (2-39).  Then 
one  has  the  polynomial 


p(w)  =  (u  +  1)  -4(u+i)  +4("+j)-K. 


(2-40) 


to  be  re-expressed  in  standard  form  as  a  polynomial  in  u.  To  do  this,  think  of  the  right  hand  side 
of  Equation  (2-40)  as  a  Taylor  expansion  of  the  polynomial  about  u  =  —  Then  express  p(u)  as 

p"(Q)  p'"(0) 

p(u)  =  p(0)  +  p'(0)u  +  ^  ■  u2  +  -  07~-u3,  and  one  has  the  canonical  polynomial  form.  One  finds  that 


3! 


ti3  -  4u/3  +  (16  -  27K)/27  =  0. 


(2-41) 


We  want  the  most  negative  root  of  the  left  member  of  Equation  (2-41),  since  that  root  will  yield  the 
least  positive  root  of  the  left  side  of  Equation  (2-39).  Using  simple  inequality  relationships,  together 
with  standard  formulas,7  one  finds  that 

«  =  1  (l  +  cos  (1  (c°s-‘  (§K  -  l)  +  2*)))  ,  (2-42) 

and  we  know  that  t2~T-  Rs2/v ,  with  s2  given  by  Equation  (2-42).  In  terms  of  K,  one  has  then 


t\—T  ^K\/G cos  f  ^Cos  1 
V  3  V  3 


=T(i-x(i+eos(KCos',( 


3tfv/6\  4tt\\ 

(2-43) 

.  )*■)) 

T-OH))) 

(2-44) 

We  give  now  a  procedure  for  obtaining  an  optimal  spiral  motion,  provided  that  v  >  R(T)/T. 
Let  (po(t),80(t))  represent  the  optimal  spiral  motion.  Then  we  want  p0(t )  =  vt,0o(t)  =  ci,0  <  t  <  tl5 
and  p0{t)  =  v{t  -  T)  +  P'T),60{t)  -  c2,<2  <  *  <  T,  where  cltc2  are  constants.  When  t  e  [<i><2],  we 
want  p0(t)  =  p{t).  From  the  relation  p2{t)  +  p2(t)02(t)  =  t;2,  we  can  secure  0(t).  In  fact,  let  f(t)  = 
dp2(t)/dt,  and  observe  that  p0(t)  =  f{t)/2po(t).  Then  we  have  p2(t)9l(t)  =  ( v 2  -  f2(t)/4p2(t)),  so  that 
0o(t)  =  ±.(4v2p2(t)  -  f2(t)) $/p2(t).  If  we  choose  the  positive  sign  for  9o(t),  we  are  assured  of  obtaining  a 
spiral  motion.  One  would,  of  course,  need  to  integrate  90{t)  numerically,  using  the  initial  condition 
0o(L)  =  ci,  in  order  to  secure  0o(t)  on  [<i,t2].  This  type  of  procedure  may  obviously  be  used  for 
any  functional  form  p(t)  whose  derivative  does  not  exceed  v  in  the  required  interval  of  t  values.  In 
particular,  when  we  develop  subsequently  the  optimal  allocation  for  preassigned  circular  disks,  one 
could  straightforwardly  obtain  the  optimal  path  in  the  same  way.  Therefore,  we  shall  not  need  to 
develop  these  paths  again  for  such  situations. 
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The  cost  functional  has  a  nonzero  integrand  on  (0, <1 )  and  on  (t2,T).  The  integral  over  (0,ti), 
using  Equation  (2-26)  and  Equation  (2-27),  is 


*-jO -£(*-=?) )’ 

-l 


dt 


**  (  v2v8t8  n2R2v6t6  irWu5t5  7r2v4ft4t4  ir Wv3R2t3 


\16a2<74  4a2<r4 


4a2<r4 


acr* 


+  W2v2t2  dt 


144a2<r4  28a2<r4 

and  the  integral  over  (t2,T)  is 


ir2v8ti  n2R2v6t \  nWvst\  ir2v4R4t\  itWv3R2t\  W2v2t\ 

4-  Trt  o  t  ~  „  o  ^  !;  o  + 


12acr2 


20a2tr4 


4acr2 


-a 


Wvt  - 


7T[t)(<  -  T)  +  R}2 
2a<r2 


(.'-W-f8!’))’ 


dt. 


(2-45) 


(2-46) 


It  is  convenient  to  make  the  change  of  variable  u  -  T -t  in  Equation  (2-46).  We  have 


Ji 


Wv(T  —  u)  — 


7 r(R—  vu)2 
2acr 2 


M)) 


du. 


The  integrand  of  the  right  hand  side  of  Equation  (2-47)  is 


,2..2,r  .a*  *Wv{T  -  u)(R-vu)2  ( r7  (R-  vuy 


;2(u)  =  W'V(T-u)2- 

+  iMzpl(R’- 

4  aJ<r4 


acr * 


( R-vu )2 


=  W2v2(T  -  u)2  -  nR2Wv{T  -  u)(R2  -  2 Rvu  +  v2u2)/aa 2 
+  7 tWv{T  -  u)(R 4  -  4 R3vu  +  6R2v2u2  -  ARv3u3  +  v4u4)/2aa2 
+  t t2R4(R4  -  4 R3vu  +  6 R2v2u2  -  4 Rv3u3  +  v4xi4)/\a2o 4 
-  t r2R2(R6  -  6R5vu  +  15flVu2  -  20R3v3u3  +  15 R2v4u4  - 
4-  v6u6)/Aa2<r4 

+  tt2(«8  -  8ft7  vu  +  28ftVu2  -  56ftVu3  +  70ftVu4  -  56 ftW 

-f  28ft2t>6u6  -  8ftu7u7  +  v8u8)/16a2cr4. 


(2-47) 


(2-48) 


Expanding  the  right  hand  side  of  Equation  (2-48),  one  finds  that 

««)  V(T  -  U)»  +  - -  UWvTS)  4 

+  "fff  lAaWvTtr1  -  ifl")  +  (fflVg1  -  .laWV1!/!  4  vT)) 
2  az(T4  2a^(T^ 

+  ^£(4<72aW(vT  +  4ft)  +  77rft4)  -  ^l(Wa3a  +  4ttR3) 

3t rW_  *2Rv7u7  ttVu8 

2a2<r4  2a2<74  16o2<t4 
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Integration  of  the  right  side  of  Equation  (2-49)  leads  to 

irR4 


7  _W2v2(T3-tl) 

Ji~  3  +  16a2 


(ir R4  -  %aWvTa2){T  -  t2)  + 


nW R4v(T  -  t2)2 
4  aa2 


",|r|’W  +  <»*»)  +  3,'flVgl- _  iIMI 


h>‘ 


12  a2c 
rr V(T  -  t2)9 


14a2<r4 


16a2  <74 


(2-50) 


144a2<r4 

We  shall  investigate  the  nature  of  both  Ji  and  J2  in  Chapter  3  when  we  compare  the  efficiency  of 
the  optimal  design  with  that  of  a  logarithmic  spiral  search. 

Let  us  consider  again  the  situation  where  effort  is  to  be  assigned  to  some  preassigned  disk  D, 
whose  probability  P(D)  is  given  a  priori.  Remember  that  Equation  (2-24)  and  Equation  (2-25)  give 
the  optimal  effort  density  assigned  to  the  disk,  depending  upon  whether  or  not  the  search  effort 
WvT  exceeds  or  does  not  exceed  the  threshold  7rc2(In2(l  -  P))/a.  We  would  like  to  know  what  the 
optimal  search  path  is  in  these  two  situations.  The  scenario  for  Equation  (2-25)  obviously  parallels 
that  for  the  unconstrained  case,  so  that  the  character  of  the  solution  just  obtained  is  preserved. 
Figure  2-7  then  applies,  with  R\{T)  replacing  R(T).  For  Equation  (2-24),  the  cumulative  allocation 
function  over  a  disk  with  center  at  the  origin  and  with  radius  r  <  R  is 


rr2  /  aW 
a  \  n 


WvT  R2  -  r2 
+ 


(2-51) 


r  R2  4  <T2 

Referring  to  our  cost  criterion  (2-26)  and  setting  c(t,p(t))  =  0,  we  find,  using  Expression  (2-51),  that 

4WvT(J2a  +  nR4  —  [(4W vTo2 a  -f  irR4)2  —  I6n  R4  cr2  aWvl}5 


p2(t)  = 


2nR2 


(2-52) 


From  Equation  (2-52),  it  follows  that 

p{t)  =  i  ^  4i4- t/jVa  *  |j  _  (1  _  \^R4<r2aWvt/{4WvT<T2a  +  */?4)2)*]  *  .  (2-53) 

Let  6(T)  =  16irft4er2aWu/(4WvTff2a+*R4)2  and  L(t)  =  b(T)t.  One  sees  that,  when  WvT  >  jr<72(ln2(l -P))/a 
(the  condition  under  which  Equation  (2-24)  applies),  4WvTa2a  >  rr/?4,  so  that  L(T)  <  1.  Obviously, 
L(0)  =  0.  Set 

ri(0  =  fl-(l  -L«»*l*  (2-54) 

so  that  ?(t)  =  (^4W R*  'j  ^  r,(t).  We  want  to  know  the  behavior  of  r,(t).  One  finds  that 

»•;(<)  =  -  a  -  L(t»*r*(i  -  L(*)r  *•  (2'55) 

From  Equation  (2-55)  rJfO)  =  +oo,  and  0  <  r[(T)  <  oo  when  4 WvTa2a  >  tR4.  Of  course,  when 
4WvTff2a  =  nR4,  rJ(T’)  =  +oo,  and  p{t)  looks  like  the  curve  displayed  in  Figure  2-2,  where  the  slope  at 
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both  t  =  0  and  t  =  T  is  +oo.  As  T  increases  such  that  iWvTcr2a  >  xRA,  the  derivative  of  p(t)  at  t  =  T 
begins  to  decrease  and,  for  a  while,  should  be  larger  than  v.  Thus,  the  character  of  our  optimal 
solution  should  be  precisely  the  same  until  such  time  T0  when  p'(T0)  =  v.  As  T  continues  to  increase, 
we  shall  show  that  p\T)  decreases,  so  that,  for  T  >  T0,  we  can  assert  that  the  radial  part  of  the 
solution  contains  only  two  parts  instead  of  three,  namely,  a  part  p  =  vt  from  the  origin  to  a  point 
(<i,p(M),  followed  by  p(t)  itself  all  the  way  from  t  =  tl  to  t  =  T.  Another  fact  will  be  of  interest:  As 
T  increases,  the  point  of  inflection  present  in  Figure  2-2  moves  progressively  to  the  right  and  finally 
disappears  leaving  p(t)  as  a  concave,  strictly  increasing  function  for  T  sufficiently  large.  We  now 
establish  these  points  mathematically. 

To  establish  the  inflectional  behavior  of  p(<),  let  us  determine  r"(<).  Using  Equation  (2-55), 
we  have 


4ri'(0  =  62m{-i[l  -  (1  -  L(t))*]-*(  1  -  L(t))-1  +  I[1  -  (1  -  !(<))*]-*(  1  -  L(<))-*},  (2-56) 

so  that 

16r«(t)  =  b\T){  1  -  (1  -  L(t))*]~i(  1  -  L(<))-^[ 2  -  3(1  -  L(t))*].  (2-57) 

From  Equation  (2-57),  r"(t)  =  0  (and  hence  p"(t)  =  0)  if  and  only  if  2-3(1  -  L(t))v  =  0.  It  follows  that 
L(hnf)  =  b(T)t, nf  =  |,  where  t,nl  means  that  time  at  which  the  curve  p(t)  has  an  inflection  point.  Note 
now  that  b(T)  is  a  strictly  decreasing  function  of  T  and  tends  to  zero  like  p  a s  T  tends  to  infinity. 
Thus  indeed  qnf  =  5/(96(T))  increases  as  T  increases  and  eventually  must  exceed  T,  so  that  the 
inflectional  character  of  p(t)  within  [0,T]  is  no  longer  present.  One  of  our  points  is  thus  established. 
We  want  next  to  show  that  p'(T )  is  a  decreasing  function  of  T,  and  we  also  want  to  find  that  value 
of  T,  say  To,  such  that  p'{T0)  =  v.  The  point  T  =  To  then  will  represent  that  value  of  T  beyond 
which  p(t)  has  only  two  parts,  one  a  nonsingular  part  and  the  other  a  singular  part  (corresponding 
to  X0  =  0).  Differentiating  Equation  (2-52)  with  respect  to  i,  we  have 


2p(t)p'(t) 


4  Wv<r2aR2 

[(4 WvTa2a  +  wR*)*  -  16irWvT<r2aR4}^  ’ 


(2-58) 


Putting  t  =  T  into  Equation  (2-58),  we  have 


mr(T) 


2  Wva2aR2 
4WvT<r2a  —  xR4 


(2-59) 


Since  p(T)  =  R ,  we  have 

~'m  -  2 Wv<r2aR 

P{  ’  \WvT<r2a  —  xR4 


(2-60) 


From  Equation  (2-60)  p\T )  is  a  strictly  decreasing  function  of  T  and  indeed  tends  to  zero  like  1/T 
as  T  tends  to  infinity.  Setting  p'(T)  =  v  in  Equation  (2-60),  we  find  that 


T0  =  R(*R3  +  2W<72a)/4Wvcr2a. 


(2-61) 
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One  last  point  to  be  made  is  that,  when  T  =  T0,tin{  <  T0.  Otherwise,  since  R  <  vT,  we  would  have 
a  logical  contradiction.  To  see  this,  refer  to  Figure  2-10.  Since  p(t)  in  this  figure  is  concave  and 
increasing,  p'(i)  >  v  for  all  t  in  (0 ,T).  Now,  since  R  <  vT0,  the  slope  of  the  line  shown  joining  the 
origin  to  ( T0,R )  is  indeed  smaller  than  v.  By  the  law  of  the  mean  for  derivatives,  there  exists  some 
t,  say  t  =  1 1,  such  that  p'(t\)  =  ^  =  kv,  where  k  <  1.  This  possibility  is,  however,  precluded  by  the 
fact  that  p'(t)  >  v  for  all  t  in  (0 ,T0).  So  Figure  2-10  presents  an  impossible  situation.  Thus,  the 
inflectional  behavior  of  the  optimal  solution  fails  to  exist  for  some  Tm{  >  T0,  where  T0,  as  before,  is 
that  value  of  t  for  which  p'(T)  =  v. 

We  next  discuss  Problem  (B):  The  target  is  stationary  in  two-space,  and  its  position  follows 
the  elliptic  Gaussian  Law  given  by  Equation  (2-10).  First  we  present  an  argument  which  will, 
simultaneously  with  the  evolution  of  the  optimal  path,  establish  the  existence  of  optimal  solutions. 
We  note  that  Hamiltonian  theory  only  affords  a  necessity  argument:  ”If  an  optimal  solution  to  the 
problem  exists,  it  is  necessary  that  certain  conditions  obtain.”  It  is  nice  to  know  that  there  indeed 
exists  at  least  one  solution  to  the  optimal  control  problem. 

P 


CONCAVE  WHEN  p'{T)  =  v 

for  Problem  (B),  instead  of  attempting  to  match  cumulative  effort  in  circular  disks,  we  now- 
attempt  to  match  it  within  ’’elliptic  disks.”  The  bounding  ellipses  become  the  level  surfaces  of  the 
density  given  by  Equation  (2-10).  Our  (local)  Lagrangian  is 

L(x,  A,  2)  =  p(x)(l  -  e~at)  -  \z,  (2-62) 

so  that  L,  =  ap(x)e~a‘  -  A.  Lt  equals  zero  when  z*  =  \  In  \  In  where  =  A/a,  as  before. 

The  optimal  solution,  analogous  to  that  for  Case  (A),  is  just  z'  =  ±  Jin  (Eiii) ] + ,  Now  p(x)  >  ^  jf 
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and  only  if  ||r||£  =  (f^-)2  +  (f^)2  <  -21n(2T<ri<r2Ai).  Thus, 


KM)  =  < 


—  ln(23r<r1o-2AJ) 


IMI2Q 


>llrll£  <  -2]n(2T£Ti<r2Ai) 

,  otherwise. 


(2-63) 


Here  ||r||^  =  rT<2r,  where  rT  =  (xi,x2)  and  Q  is  the  diagonal  matrix  Q  =  diag(l/crf,  The  total 

search  cost  is  again  WvT,  where  v  is  the  speed  in  the  (xi,x2)-plane.  Integrating  /^(x)  over  x  space, 
one  finds  that 

WvT  =  ^^■ln2(2x<r1<72A1),  (2-64) 

a 

so  that  Equation  (2-63)  becomes 


[l 

/aiwry 

*  IM|2q' 

/:,(*)  =  - 

1  a 

1 

y  X(Ti<T2  / 

2 

l  0 

.otherwise. 


(2-65) 


Integrating  the  right  hand  side  of  Equation  (2-65)  over  the  elliptic  disk  D:  {x  :  (f^)2  4-  (f^)  <  r2  < 
J?2(T)},  where  R2(T)  =  2  ( aJ%f£, )  ^ 1  one  fin<^s  that  the  cumulative  effort  £(||r||<j)  is 


E(Mq) 


2  a 


Our  cost  functional  is  now 


(2-66) 


J  =  C(r(t))=  fT(c(t,r{t)))2dt 
Jo 

=  [  (Wvt  —  £(||r(<)||<j))2  dt,  (2-67) 

Jo 

with  £■  (||r(T)||Q)  afforded  by  Equation  (2-66).  The  reader  may  check  the  fact  that  E  (|[j*(T,)[|q)  = 
E(R(T))  =  WvT,  the  total  search  effort,  provided  we  require  ||r(T)||Q  =  R(T).  The  searcher  motion  is, 
of  course,  given  by  r(<)  =  v,  where  the  usual  (Euclidean)  norm  of  v  ,  namely,  ||v||,  is  just  the  speed 
t.  Furthermore,  we  want  r(0)  =0,  and  we  desire  that  ||r(t)||Q  be  monotone  increasing;  thus 


r  Qv  >  0. 


(2-68) 


The  problem  which  we  want  to  address  is  that  of  determining  whether  or  not  there  exists  a 
solution  to  the  problem 

minC(r(<)),  (2-69) 

where  r(0)  =0,  r(<)  is  piecewise  continuous  on  {0, T],  [|r(£)|j  =  v  whenever  r(t)  is  defined,  and  ||r(r)||o  = 
R(T),  where  Q  is  any  positive  definite  symmetric  matrix.  Note,  of  course,  that  the  standard  Eu¬ 
clidean  norm  is  used  on  the  derivative  of  r (i).  We  cannot  state  categorically  that  such  a  problem 
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will  have  a  solution.  On  the  other  hand,  if,  instead  of  requiring  that  ||r(<)||  =  v  except  for  a  finite 
number  of  values  of  t,  we  demand  only  that  ||r(<)||  <  *>,  standard  theory8  indicates  that  a  solution 
can  indeed  be  found. 

Let  Xo  be  the  set  of  functions  r(i)  in  the  domain  of  our  original  problem  and  X\  be  the  set 
of  r (<)  for  the  new  problem.  Clearly  ^  C  X).  We  claim  that  X0  is  dense  in  Xx  in  the  sup  norm 
topology.  That  is,  given  any  f(t)  in  X i  and  t  >  0  arbitrary,  but  positive,  there  exists  f((<)  in  Xo  such 
that 

sup  ||f(t)-fc(*)ll<‘- 
0  <t<T 

The  situation  is  illustrated  in  Figure  2-11  , which  shows  a  "zigzag  function  £<(0”  approximating  f (<). 
Here  f<(t)  =  f(t)  for  0  <  t  <  ti,  and  f{(<)  #  f (t)  on  <i  <  t  <  T.  Clearly,  the  amplitude  and  frequency  of 
the  zigzag  can  be  so  controlled  that  the  searcher  speed  is  precisely  t>  at  all  points  of  the  path  r  =  f e(t) 
and  that  ||f4(t)  -f(t)||  <  t  for  all  t  in  (0,7''].  Furthermore,  the  following  assertion  is  now  obvious: 

inf  J  =  min  J. 
r(t)eXo  r(*)€X, 

We  don’t  know  a  priori  that  J  attains  its  infimum  on  but  we  can  examine  the  minimizing 
solution  in  Xi  and  (if  it  is  not  in  X0)  approximate  it  by  a  point  of  X0  (i.e.,  a  function  of  t)  as  closely 
as  we  like. 


X2 


FIGURE  2-11.  ZIGZAG  FUNCTION  APPROXIMATING  A 
GIVEN  FUNCTION 
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We  asserted  previously  that  r(t)  is  piecewise  continuous  on  [0,  T’J,  i.e.,  its  components  i,-(*), 1  = 
1,2,  are  piecewise  continuous  on  that  interval.  It  can  then  be  shown  that  p(t)  =  ||r(t)||<3  is  a  Lipschitz 
function  and  that,  therefore,  it  is  certainly  absolutely  continuous.  9  We  have,  for  any  tltt2  belonging 
to  [0,  Tj, 


It  remains  to  find 


|p(*i)  —  /K*2)l  =l  llr(*i)ll<?  —  |lr(<2)|lo  | 

<  lk(h)  -  r(<2)||<3 

<\  f 3  \\m\Qdt\ 

Jti 

<  |  [  max  ||r(<)ll<3  dt\ 


(2-70) 


max  ||r(<)ll<?.  (2-71) 

IW«)II<* 

i.e.,  the  maximum  of  the  Q-norm  of  the  velocity  vector  when  its  ordinary  Euclidean  norm  does 
not  exceed  v.  In  other  words,  we  want  to  find  a  vector  x0  whose  usual  norm  does  not  exceed  v 
and  which  maximizes  the  expression  x'Qx.  Now,  since  Q  is  symmetric,  there  exists  an  orthogonal 
matrix  P10  such  that  Q  =  P'DP,  where  D  is  a  diagonal  matrix  and  the  entries  on  the  diagonal,  say 
l/<7jf  and  l/<rf,  are  the  eigenvalues  of  matrix  Q.  Assume  that  <r2  <  ci-  The  form  x'Qx  is  y'Dy  in  y 
coordinates,  where  y  =  Px.  Since  P  is  orthogonal,  it  is  well-known  10  that  ||y||2  =  ||x||2  in  the  usual 
sense,  so  that  ||y||2  <  v.  Now 


y'  Dy  =  y\/o\  +  y2/ff2 

<  iy\  +  y\)l°\ 

<  v2  jcr\. 


Clearly,  the  vector  y0  =  (0,v)T  in  the  y-coordinate  frame  is  a  candidate  y  which  maximizes  our 
quadratic  form.  This  vector  is  an  eigenvector  of  length  t;  corresponding  to  1  /<r\  for  the  diagonal 
matrix  D.  Therefore,  x0  =  P'yo  is  likewise  an  eigenvector  of  Q  for  the  eigenvalue  1  /<r\.  It  follows 
that,  whenever  r(<)  exists  (it  is  indeed  piecewise  continuous),  its  Q-norm  is  largest  when  r(l)  is 
an  eigenvector  of  Q  corresponding  to  1  /a\.  This  would  occur,  in  particular,  when  r(<)  is  on  the 
semiminor  axis  of  the  ellipse  E  =  { y  :  y2/<r2  +  y\l°\  —  1} ,  r(t)  points  along  the  semiminor  axis  away 
from  the  origin,  and  ||r(i)ll2  =  v.  Thus,  from  Inequality  (2-70),  we  have 

\p{U)~  p(t2)\  <  v|<2  -  ti\/<r2,  (2-72) 

so  that  p(t)  is  Lipschitz  continuous  on  [0,  T\  with  Lipschitz  constant  v/a2.  By  Lipschitz  continuity, 
p{t)  exists  almost  everywhere  on  [0,  T]  and  is  Lebesgue-summable,  the  integral  of  p(t)  being  p{t)  itself. 
From  Inequality  (2-72)  we  discover  that,  wherever  p(t)  exists ,  it  is  bounded  in  absolute  value  by 
v/p2,  i.e.,  |p(<)|  <  d/p?  almost  everywhere  on  [0,7^. 
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Suppose  that  the  c  =  0  path  qualitatively  appears  as  shown  in  Figure  2-12.  We  shall  show 
later  that,  for  our  effort  function  given  by  Equation  (2-66),  this  assertion  is  valid.  Then,  in  analogy 
with  Case  (A),  consider  the  path,  as  shown,  which  moves  out  from  the  origin  along  a  semiminor  axis 
at  speed  v(p  =  v/a2)  until  c  =  0,  whereupon  t  =  U,  then  follows  the  c  =  0  curve  at  a  speed  less  than  v 
while  still  moving  along  the  semiminor  axis  until  a  certain  time  t2,  whence  it  switches  back  to  speed 
v,  continuing  to  move  along  the  semiminor  axis  up  to  the  terminal  time  T,  at  which  ||r(T)||Q  =  R(T). 
The  switching  time  t2  is  determined,  as  shown  in  Figure  2-12,  by  the  point  of  intersection  of  the  line 
p  =  £(t  -T)  +  R(T)  and  the  curve  c  =  0 (p  =  p(t)).  Note  that,  at  the  moment ,  we  are  only  considering 
a  path  which  coincides  with  the  semiminor  axis  of  the  ellipse  ||r||Q  =  R(T).  Let  r.  be  the  above  path, 
with  corresponding  Q-norm  p.,  and  compare  it  with  any  other  admissible  path  r  (with  Q-norm  p). 

From  our  above  analysis,  the  maximum  value  of  p(t),  when  it  exists,  is  t >/<r2,  so  that,  almost 
everywhere  on  both  [0,<i]  and  [t2lT\,  p(t)  <  v/o2  =  p,(t).  Consider  the  interval  [0,*i].  Let  f(t)  = 
p(t)  -  p,(t)  there,  and  we  see  that  f(t)  <  0  almost  everywhere.  Thus,  for  any  <6[0,<i], 


i: 


f(r )  dr  <  0, 


so  that 


m-m<o. 


But  /( 0)  =  0,  since  both  paths  must  emanate  from  the  origin.  It  follows  that  /(<)  <  0  on  [0,<i], 
i.e.,  that  p(i)  <  p.(t)  there.  On  the  other  hand,  we  also  know  that,  for  t  in  [t2,T],f(t)  <  0  almost 
everywhere,  so  that 


l 


f(r )  dr  <  0. 


Thus, 


f(T)  -  f(t)  <  0. 

Now  f(T)  =  0,  since  both  paths  must  satisfy  the  condition  ||r(T)||q  =  R(T).  So  f(t)  >  0  on  [<2,T], 
implying  that  p(t)  <  p,(t)  there  ience,  appealing  to  Figure  2-12  again  and  to  the  nature  of  c(t,r(<)), 
we  see,  for  example,  that 

c(t,  ?(<))  >  c(t,  r.(t))  >  0,  t  €  [0, «!].  (2-73) 

The  first  of  inequalities  (2-73)  follows  because  p(t)  <  p.{t)  in  [0,<i],  so  that  E(p(t))  <  E(p.(t))  there. 
Thus  c(t,r(())  =  Wvt~  E(p(t))  >  Wvt  —  E(p.(t))  =  c(t,rm(t)).  The  second  inequality,  namely,  c(t,r.(t))  >  0, 
is  vahd  because  r.(<))  lies  in  the  region  c  >  0  of  Figure  2-12.  Also,  by  definition  of  r.(r), 


|c(<, r(<))|  >  |c(t, r.(<))|  =  0,<  €  [t\,t2]- 
Finally,  by  reasoning  paralleling  that  for  Inequalities  (2-73), 

c(f,r(0)  <c(Lr.(<))  <  0 ,«€  [t2,T\. 


(2-74) 


(2-75) 
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It  follows  that  C(t,r(t))  >  C(t,r.(t)),  i.e.,  that  r,  is  optimal  in  the  class  X\. 

As  illustrated  in  Figure  2-11,  we  can  approximate  r,  by  a  path  in  X0  (for  which  the  speed  is 
always  v )  as  closely  as  we  like.  (As  we  previously  stated,  a  zigzag  function  will  do.)  Now  if,  instead 
of  a  zigzag  path,  we  can  arrange  a  spiral  path  in  X0  which  connects  (*i,r.(ti))  to  (<2 , r. (<2))  while  also 
maintaining  c  =  0,  we  will  have  an  optimal  path  in  Xo .  In  other  words,  we  can  certainly  obtain  an 
e-approximate  solution  in  Xo  to  an  optimal  solution  in  Xi,  but  there  is  some  doubt  that  a  minimizer 
exists  in  X0  for  any  choice  of  the  pair  (<71,02).  However,  it  is  easily  shown  that,  when  0^  =  <r2,  a  spiral 
solution  which  meets  the  requirements  can  be  generated  if  only  R  <  vT.  One  might  conjecture  then 
that,  when  o^R  <  vT,  an  elliptic  spiral  meeting  the  conditions  could  likewise  be  obtained  (  the  speed 
being  larger  that  that  required  to  proceed  from  the  origin  along  the  semimajor  axis  of  the  bounding 
ellipse,  arriving  at  the  boundary  at  time  T).  In  any  event,  for  every  pair  (<ti,<72),0  <  <r2  <  0^,  there 
hopefully  exists  vo(<ri,02,  T)  such  that,  when  v  >  vo(0’i,a2,7'),  spiral  motion  is  possible. 


Pit) 


FIGURE  2-12.  OPTIMAL  SEARCHER  VERSUS  NONOPTIMAL 
SEARCHER  MOTION  IN  ELLIPTIC  CASE 

We  now  show  that  the  c  =  0  curve  in  Figure  2-12  has  the  qualitative  character  depicted  there. 
In  fact,  it  is  possible  to  prove  that  p(t)  is  exactly  the  same  function  as  that  obtained  for  Case  (A), 
i.e.,  that  Equation  (2-28)  holds,  with  p(t)  =  ||r(<)J|g.  Indeed,  from  Equation  (2-66)  and  from  the  fact 
that  p(0)  =  0 

*,F(<»  =  <Vv<  -  =  0 
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if  and  only  if 

p2(t)  =  R7(T)[1  -  (1  -  4aWvt/ir<rl<r2R4(T ))*].  (2-76) 

Now  R?(T)  =  2(aWvT/ir<Tia2)l,  so  that,  from  Equation  (2-76),  p7(t)  =  R7(T)[  1  -  (1  -  t/T)i ),  and  we 
again  have  Equation  (2-27),  but  this  time  with  an  interpretation  for  p(l)  consistent  with  an  elliptic 
norm. 

Just  as  in  the  circular  case,  it  is  possible  to  design  a  procedure  for  generating  an  optimal 
spiral  search  pattern,  provided  i>0(<7i,<t2,T)  exists  for  any  («71,a2)-  The  procedure  (algorithm)  which 
we  now  present  should  allow  one  in  practice  to  devise  the  path  if  it  is  possible  to  do  so.  Let 


*?(<) 


+ 


!W  = 


=r(t)  =  /(<), 


(2-77) 


and  let  us  require  that 


*?(0  +  *j(0  =  *2- 


(2-78) 


Making  the  transformation  yi(t)  =  *i(t)/<rx,y2(t)  =  x2(*)/c2  in  Equation  (2-77)  and  Equation  (2-78), 
we  find  an  equivalent  pair  of  relations,  namely, 


yi(<)  +  »|(0  =  /(<) 


(2-79) 


and 

*i>?(0  +  *liS(0  =  »a-  (2-80) 

Let  r7(()  =  y?(<)  +  y|(<)  and  yx (<)  =  ry(<)cos0y(f),y2(t)  =  r2(<)sin0y(<).  Then  one  sees  that 


(T7 (r7(t)  cos7  6y(t)  -  ry(t)ry(t)(sin20y(t))0y(t) 

+ry(0(sin2  6y  (t))07(t))  +  a\  (r7(t )  sin2  0y{t) 

+ry(OMO(sin20y(t))Mt)  +  rJ(t)(cos2  0,(O)*?(<))  =  v7.  (2-81) 


Now  we  use  the  fact  that  fy(t)  =  /(t)/2ry(f),  substitute  this  relation  into  Equation  (2-81),  and  solve 
the  quadratic  equation  in  0y(t)  to  obtain 


where 


and 


g  M  _  /(<)(^?-gI)Bin(2gy(0)±2/l(Q) 
V(  }  4 /(!)/»(!) 


/i(0s(4»2/(0/jW-^!/,(0)1 


/2(i)  =  ff2sin20y(t)  +  erf  cos2  0y(t). 


(2-82) 


To  insure  that  0y(t)  is  positive,  let  us  choose  the  positive  sign  in  Equation  (2-82).  Using  the  initial 
condition  0y(<i)  =  »/2,  we  could  then  solve  Equation  (2-82)  numerically  to  obtain  0y(«),  <,  <  t  <  t7.  Now 
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*i(0  =  =  ffir„(t)co6  0j,(*),x2(t)  =  o2y2{t)  =  <T2ry(<)sm5y(t),  and  r|(«)  =  r*(t)(<7 \  cos2  0s(*)+<x2  sin2  0y(<))- 

Also,  0X(*)  =  tan-1(ar2(i)/xi(t))  =  tan_1(<r2  tan0y(t)/<r1).  Since  the  inverse  tangent  and  tangent  functions 
are  both  strictly  increasing  functions  of  their  arguments,  0x(t)  is  an  increasing  function  of  t  whenever 
0y(<)  is.  The  pair  (rx(<),0x(<))  thus  obtained  would  represent  an  optimal  solution  on  [<i,t2],  given 
rx(ii)  =  v<i,0r(ti)  =  jt/2  as  initial  information.  The  solution  found  thus  far  would,  however,  not  be 
satisfactory,  since  we  require  that,  at  time  t2,(x i(t2),x2(<2))  is  to  be  on  the  semiminor  axis  of  the 
bounding  ellipse.  To  rectify  this  situation,  we  start  from  the  point  rx(t2 )  =  a2R{T)  -  v(T  -t2),0z(t2)  = 
?r/2  on  the  semiminor  axis  and  allow  9X  to  increase  as  time  decreases  away  from  t2.  To  implement 
such  a  procedure,  again  choose  the  positive  sign  in  Equation  (2-82).  Now  let  (r^(t),0/(<))  be  the 
forward  solution  obtained  by  marching  from  t  =  ti  forward  to  t  =  t2,  and  let  (r£(<),0x(O)  be  the 
backward  solution  obtained  by  receding  from  t  =  t2  to  t  =  ti.  Then,  if  the  forward  and  backward 
solutions  intersect  at  some  point  i  =  tiNT,  an  optimal  solution  on  [<i,<2]  would  be  (r"(<),0°(*)),  where 
(r2(0.*2(0)  =  (ri(0.^(<)).*i  <  <  <  iiNT,  and  r”(t),0‘(<))  =  (r£(t),0£(<)).</WT  <  t  <  h-  It  is  possible  to 
show  that  they  will  intersect.  To  see  this,  examine  first  the  circular  situation.  Refer  to  Figure 
2-13.  If  we  consider  the  total  angular  increment  from  time  t\  to  time  /2,  divide  the  increment  in 
half,  and  then  reflect  the  portion  of  the  curve  from  point  (r(th),0(tfc))  to  (r(<2),0(<2))  about  the  line 
shown,  we  have  an  equivalent  path  with  respect  to  allocation  considerations.  In  so  doing,  the  point 
(r(<2),0(<2))  gets  mapped  to  (r(<2),  t/2).  We  can  do  this  because,  in  the  circular  case,  the  allocation 
density  is  radially  symmetric.  Let  and  e2  be  small  and  positive.  Imagine  that  one  were  to  change 
the  circle  of  radius  R(T)  into  an  ellipse  whose  semiminor  axis  is  aligned  with  the  x2  direction  and 
has  length  (1  -  c i){R(T )  and  whose  semimajor  axis  is  aligned  with  the  x,  direction  and  has  length 
(1  +  (2)R(T).  Then,  by  continuity,  the  modified  optimal  path  Pc  (e  for  circular)  would  be  converted 
into  an  optimal  path  Pe.  That  path  Pe  must  start  on  the  semiminor  axis,  namely,  the  x2  -axis,  and 
end  on  the  semiminor  axis,  namely,  the  x2  -axis  again.  The  backward  trajectory  portion  of  Pe  will 
simply  be  a  modification  of  the  dashed  part  of  Pc  leading  from  (rl(t2),v/2)  back  to  (r(<A),0(t*)),  and 
the  forward  part  of  Pe  will  correspond  to  the  forward  portion  of  Pc.  Therefore,  such  an  intersection 
point  is  expected  to  exist,  in  general;  since,  as  the  ellipticity  is  progressively  changed,  maintaining 
the  x2  -axis  as  the  sem.uiinor  axis,  the  optimal  path  must  change  in  such  fashion  that  one  starts  on 
the  x2  -axis  and  terminates  on  the  same  axis. 

The  argument  here  is  not  a  proof  of  existence  of  an  optimal  solution  in  X0.  It  should  be 
taken  rather  as  a  plausibility  argument.  We  are  using  what  mathematicians  call  a  homotopy  in 
order  to  derive  optimal  solutions  for  new  situations  from  those  obtained  under  old  conditions  which 
are  close  to  the  ne  v  ones.  One  may  use  such  a  concept  to  obtain  numerical  solutions  to  Equation 
(2-82),  employing  the  so-called  Picard  method.11  As  input  to  the  method,  one  might  start  with  0y(<) 
obtained  for  the  circular  case. 

Parallel  to  the  circular  case,  one  may  assume  that  all  effort *is  to  be  devoted  to  a  preassigned 
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FIGURE  2-13.  ALTERATION  OF  PATH  IN  CIRCULAR  CASE 
TO  MEET  PATH  CONSTRAINT 

elliptic  disk  Dt.  Suppose  that  we  represent  this  disk  by  De  =  {(xi,x2)  :  (i\/crl)2+  {x2/a2)2  <  R 2).  Then 
we  have 


P(D-}  =  IJ^b;'xp  H  (?f +  ® )  dx‘ dxi- 


(2-83) 


Letting  yi  =  x1/(rl,y2  =  x2/<?2,  Equation  (2-83)  becomes 


p{°e )  =  JJ  exp(-^(y?  4-  y\))dy\dy2, 


(2-84) 


where  D  is  a  circular  disk  of  radius  R.  If  we  now  change  to  polar  coordinates,  i.e.,  let  yx  -  rcos0,y2  = 
r  sin  6,  Equation  (2-84)  becomes 

P{De)  =  (I  —  exp(-/22/2)).  (2-85) 

From  Equation  (2-85)  we  have 


R=  [-21n(I-P(D<))]i1 
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so  that,  given  P(D«),  one  knows  R  and  hence  the  disk  De.  Parallel  to  Equation  (2-20),  one  has  the 
conditional  Gaussian  distribution 


p„.(x)  =  /  (-5  (I  +  If)) W. 

l  0 


)  ,  x  6  De 

.x  £  De. 


(2-86) 


We  use  Equation  (2-21)  with  po(x)  replaced  by  pd.(x)  in  Equation  (2-86).  As  before,  let  A!  =  A/a. 
Then  pD.(x)  >  A!  if  and  only  if  ||x||£  <  -21n(2xcr1(T2A1P).  So 


-  \  -\n(2v(Tl(T7XiP(De))  - 


n,oM)  =  •  ["  -  INI h  < 

•0  .IWIq  > 

where  R *  =  R?(P, Ax,^,^)  s  min(J?2,-2b(2T<r1<r2A,P)).  Now  define  A2  =  A ,P(De),  and  compute  the 
cost  over  the  elliptic  disk.  We  have 


ClKoM))  = 


Per i  a?  R}  1  i  _  i 

— —  ln(2x<Ti<72A2) — -  ,0<A2<~ 

“  L  4  l*(T\  C 


ir<ri<r2 


fin(2jr<ri<r2A2)]2 


1-P 

’  2x<Ticr2 


(2-87) 


<  A2  < 


2k<T\<T2 


The  first  line  of  Equation  (2-87)  corresponds  to  the  case  where  P2(P,A,,<71,<t2)  =  R2.  Then  /j  D  is 
nonzero  over  the  entire  elliptic  disk  ||x||Q  <  R.  Equating  the  right  hand  side  of  line  1  of  Equation 
(2-87)  to  A,  =  WvT ,  one  finds  that  A2  =  ^exp  (-£  -  Inserting  this  relation  for  A2  into 

0  <  A2  <  (1  _  P)/2x<ti<t2,  we  find  that  A,  >  x(ln(l  -  P))V,<r2/a.  Also,  substituting  the  A2  value  obtained 
above  into  the  allocation  density,  we  have 


n.oM)  =  i  ( 


attA>T  fl2  -  2||x|l^ 
TtR2a\<T2  4 


,1|x||q  <  R,WvT  >  x<Ti(72(In2(l  -  P))/a. 


(2-88) 


In  addition,  note  that  a  searcher  can  reach  the  boundary  of  an  elliptic  disk  after  T  units  of  time  if 
and  only  if  he  can  do  so  along  its  semiminor  axis.  Thus  we  should  require  that  Ra2  <  vT,  i.e.,  that 
R  <  vT/<t2. 

Now  consider  the  second  line  of  Equation  (2-87).  Then  0  <  A,  <  jr()n2(]  -  P))ffi<r2/a,  corre¬ 
sponding  to  2J~^a  <  A2  <  Also,  A2  =  2di<73  exp  -  •  Inserting  this  equality  into  the 

allocation  density,  we  have,  when  WvT  <  *(ln2(l  -  P))«r,<r2/a, 


/a,d,(x)=  2  "  "Q  -  1  - 

lo  ,R\  <  ||x|l^  <  P2, 


(2-89) 


where  fi2  —  -21n(2x<r1<r2A2)  =  2(aWvT/x(ri<T2)a .  Equation  (2-89)  has  the  same  form  as  Equation 
(2-65),  the  difference  between  the  two  being  that  Equation  (2-89)  presumes  a  distribution  over  a 
preassigned  elliptic  disk  and  that  Equation  (2-65)  presumes  one  over  the  entire  x  plane.  When 
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the  effort  to  be  expended  does  not  exceed  *-(ln2(l  -  P))tr\<j2/a,  the  situation  is  parallel  to  that  of 
unrestricted  disk  size.  The  theory  just  formulated  for  Case  (B)  then  dictates  the  optimal  path. 
The  optimal  path  process  for  Equation  (2-88)  remains  to  be  generated.  We  shall  next  develop  this 
process. 

The  cumulative  allocation  function  (cumulative  effort)  over  a  subdisk  ||xj|<}  <  r  of  the  elliptic 
disk  ||x||q  <  R  when  WvT  >  ir0-i<T2(ln2(l  -  P))/a  is  found  by  integration  of  Equation  (2-88)  over  that 
elliptic  subdisk.  We  find  that 


E(Mq)  =  E(p)  = 


/  aWvT  R2  -  p2\ 
\irR2<rl<T2  +  4  / 


,p  <  R. 


(2-90) 


Referring  to  our  cost  criterion  (2-67)  and  setting  c(t,r(t))  =  0,  we  find,  using  Equation  (2-90)  and 
the  fact  that  p(0)  =  0,  that 


p\t)  = 


AaWvT  -f  iraiazR*  —  [(4aWvT  +  ncri^R4)2  —  l&iracriaiR^Wvt} 5 

2*  (TitT^R2 


From  Equation  (2-91)  it  follows  that 


(2-91) 


_  1_  / 4aiy«T  +  ■k<t1(T2R?\^  /  lGnaci^R^Wvt  ’ 

P[  ~  R  V  2jror1op2  )  ~  V1  "  (4aWvT  +  ira^R4)2 ) 


(2-92) 


Let  6(T)  =  \S*a<TXo-,R?Wvl{\aWvT +ir<71<r2/i4)2  and  L(t)  =  b(T)t.  One  sees  that  WvT  >  waxa2(^2(l-  P))/a 
if  and  only  if  AaWvT  >  x*i<r2&,  so  that  L(T)  <  1.  Obviously,  L{ 0)  =  0.  Set 


r,(0=  [l-(l-L(0)^,  (2-«3) 

the  same  form  as  Equation  (2-54).  Then  p(t)  =  -fa  ri(0-  Just  as  in  the  circular 

case,  we  can  assert  that  ^(0)  =  +oo  and  that  0  <  r[(T)  <  oo  when  our  new  condition  WvT  > 
jro-1<T2(ln2(l  -  p))/a  applies.  Of  course,  when  4 aWvT  =  * o-i^P^.rjlT)  =  +oo,  and  p(t)  looks  like  the 
curve  displayed  with  hatch  marks  in  Figure  2-12.  As  T  increases  such  that  4 aWvT  >  nai(T2R4,  the 
derivative  of  p(t )  at  t  =  T  begins  to  decrease  and,  for  a  while,  is  larger  than  v.  From  the  right  side  of 
Equation  (2-55),  which  is  also  symbolically  the  derivative  of  the  right  side  of  Equation  (2-93),  we 
see  that  r[(T)  tends  to  zero  as  T  tends  to  infinity.  We  shall  show  that  p'(T)  is  a  strictly  decreasing 
function,  so  that,  at  precisely  one  value  of  T,  say  T  =  T0,p'(T0)  =  v/<r2.  For  T  >  T0,  the  line  segment 
portion  of  p*(t)  in  Figure  2-12  disappears,  and  it  becomes  possible  to  follow  the  c  =  0  curve  all  the 
way  from  t  =  <i  to  t  =  T.  The  inflectional  nature  of  rx(t)  is  precisely  the  same  as  in  the  circular  case, 
and  the  analysis  presented  there  can  be  repeated  verbatim  here.  Indeed  tmt  =  5/(96(T))  again,  with 
6(7^  a  strictly  decreasing  function  of  T,  behaving  like  1  /T2  as  T  tends  to  infinity.  Thus,  as  in  the 
circular  case,  the  inflectional  character  of  p(t)  is  no  longer  present  when  T  is  sufficiently  large. 
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We  would  now  like  to  show,  as  in  the  circular  case,  that  p'(T)  is  a  strictly  decreasing  function 
of  T,  and  we  want  to  ascertain  T0,  where  p'(T0)  =  «/<r2.  From  Equation  (2-91),  we  have 


p(T)p\T)  = 


2aWvR2 

AaWvT  —  mriffiFP 


(2-94) 


Since  p(T)  =  R(T),  one  finds  that 


P'(T)  = 


2  aWvR 

AaWvT  —  TOiO-iP* 


(2-95) 


Thus  p'{T),  as  in  the  circular  case,  is  a  decreasing  function  of  T,  tending  to  zero  like  l/T.  Setting 
p'(T)  =  vja-x  in  Equation  (2-95),  we  have 


+  2aW) 

0  ~  AaWv 


(2-96) 


Also,  similar  to  the  argument  used  for  Figure  2-10,  we  can  claim  that  T,n f  >  T0,  so  that, 
whenever  T  >  Tinf,p(<)  would  appear  as  shown  in  Figure  2-10,  where,  of  course,  for  the  elliptic  case, 

m  =  IWOIIq- 

Parallel  to  Case  (A),  we  would  next  like  to  consider  the  order  of  the  error  when  one  substitutes 
p.(t)  into  Equation  (2-67).  (We  shall  only  investigate  the  error  for  the  unconditional  elliptic  Gaussian 
distribution,  i.e.,  the  one  corresponding  to  a  moving  boundary  ellipse.)  On  0  <  t  <  t\,p.(t)  =  vt/a7; 
using  Equation  (2-66)  and  Equation  (2-67),  we  have  the  following  parallel  to  J i  for  the  circular  case 
(given  in  Equation  (2-45)): 

_  r2<r2v8tf  7T2R2v6tj<T2  w<riWv5tf 

144a2cr§  28aV^  I2aa% 

JOxWv3R2t\  W2v2t3 
20a2 a2  4a<r2  3 


For  J2,  one  has 


-n 


Wvt  - 


ir<T\<Ti  ( R  +  v(t  -  T)!<rj)‘ 


R2  - 


(R  +  v(t  -  T)/a2y 


Wv(T  -u)-  ™1<T2  (f?2~  vuJ°-2i-  ^_(g~^)2  \  du 


(2-98) 


upon  making  the  change  of  variable  u=T-t.  Expanding  the  integrand  of  the  right  member  of 


2-30 


NAVSWC  TR  87-306 


Equation  (2-98)  and  integrating  the  result  term  by  term,  we  find  that 

W2  v2  (X3  —  <2)  *o\Oi  R?  (*010210*  —  SaWvT)  (X  —  (2) 

h  = - ? - +  — ~ ( - ik - 


+ 


+ 


WvR*  (X  -  <2)2  R?Wv3  (X  -  <2 )4 


2iT2 


td  (7n<Tl<r2R4  +  4aW  (4 fig;  +  vX))  (X  -  t2)5 


40ao-2 


+ 


+ 


t;5  (4*<7ifl3  + Wa)(X-<2)6 

12a<x4 

3T<Tii?2t/6  (X  —  t2f  *<J\Rv2  (X  —  t2)8 


14  a<r| 


16a^| 


ircriu8  (X  —  <2)9 


144a<T2 


} 


(2-99) 


The  values  for  ti  and  t2  are  again  obtainable  through  Equation  (2-43)  and  Equation  (2-44), 
but  with  a  different  definition  for  R(T).  That  is,  Equation  (2-43)  and  Equation  (2-44)  yield  and  t2, 
respectively,  when  K  =  ^  and  R2(T)  is  given  by  2  a  instead  of  2 a2  (°-y~r-)  *  =  2<r  (°~~-)  *  •  We 

shall  investigate  the  nature  of  Ji  and  J2  as  X  tends  to  infinity  in  Chapter  3,  which  we  turn  to  next. 
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CHAPTER  3 

COMPARISON  OF  THE  EFFECTIVENESS  OF  OPTIMAL  CONTROL 
STRATEGIES  WITH  LOGARITHMIC  SPIRAL  DESIGNS 

Up  to  this  point,  we  have  developed  the  optimal  allocations  of  effort  for  stationary  targets 
whose  distributions  are  known  to  be  Gaussian,  and  we  have  devised  procedures  for  generating 
optimal  search  strategies  when  these  are  expected  to  exist.  Optimal  search  plans,  even  if  they  turn 
out  not  to  be  practically  useful,  serve  as  a  basis  of  comparison;  after  all,  they  represent  the  best 
that  one  can  do  under  the  stated  conditions  of  the  problem.  The  purpose  of  this  chapter  is  (1)  to 
develop  logarithmic  spiral  solutions  to  our  problem  for  both  the  circular  and  elliptic  Gaussian  cases 
and  (2)  to  compare  the  asymptotic  properties  of  the  optimal  solutions  developed  in  Chapter  2  with 
those  for  the  spirals.  We  shall  discover  that  the  optimal  solutions  yield  a  mean-square  error  of  the 
order  of  T~$  as  T  tends  to  infinity,  v  being  held  constant,  and  that  the  logarithmic  spiral  solutions 
behave  like  T3  in  this  regard.  On  the  other  hand,  we  shall  find  that,  if  we  properly  choose  v  as  a 
function  of  T,  we  can  drive  both  the  cost  for  the  spiral  and  the  optimal  solution  toward  the  same 
asymptotic  behavior,  thus  making  the  logarithmic  spiral  computationally  attractive. 

Let  us  begin  by  studying  the  properties  of  a  conventional  logarithmic  spiral,  which  we  might 
like  to  adapt  to  the  situation  where  the  target  is  stationary  and  its  location  follows  a  circular 
Gaussian  law.  The  general  form  of  the  logarithmic  spiral  is  p  =  cue'3*,  where  cj  and  c2  are  given 
positive  quantities,  which  may  depend  upon  the  budget  time  T,  the  speed  v  of  the  searcher,  the 
searcher’s  detection  width  W ,  and,  finally,  the  variance  o-2  of  the  target’s  location.  It  is  important  to 
note  that  <u  and  c2  are  to  be  constants  once  such  parameters  are  adjusted.  The  problem  of  interest 
thus  becomes:  Find  p  and  6  as  functions  of  time  t,  where  0  <  t  <  T,  such  that  p  =  cue'3*,  the  speed 
along  the  spiral  is  v ,  ?nd  the  searcher  meets  the  boundary  of  a  known  circular  disk  p  <  R(T)  at 
time  t  =  T.  Then  compute  the  cost  function,  given  by  Equation  (2-26),  where  p  =  p(t)  is  the  radial 
component  of  the  logarithmic  spiral. 

Our  target  distribution  for  the  circular  Gaussian  case  is  given  by  Equation  (2-9).  w'here 
r2  =  x2  +  of  course.  Since  the  target  is  most  likely  to  be  at  or  close  to  the  origin  of  the  (au,x2) 
plane,  we  agree  to  start  the  spiral  at  the  origin  (or  perhaps  at  i  distance  away  from  it)  and  spiral 
out  to  the  boundary  of  the  disk  p  <  R(T ),  on  which  the  optimal  effort  happens  to  be  supported. 
The  cumulative  effort  to  be  approximated  is  given,  in  fact,  by  Equation  (2-27).  Note  that,  w'hen 
6  =  -oo,  p-  0,  so  that  one  has  the  initial  information:  At  time  t  =  0,  p  =  0  and  0  =  -oo.  We  shall  see 
that,  when  speed  is  the  sole  requirement  to  be  matched,  one  can  design  a  spiral  which  is  traversable 
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at  that  speed  from  the  origin  to  the  boundary  of  disk  p  <  i£(T),  reaching  the  boundary  at  time  T, 
provided,  of  course,  that  R(T)  <  vT. 

Now  the  length  of  a  curve  p  =  f(9)  between  two  points  (pi,0i)  and  (p2,0 2),  where  p  and  6  are 
polar  coordinates,  is  given  by12 


L  =  ['  (f2{6)  +  f\6)y  d9, 

j  6  D 


so  that,  when  p=  cieC30,ci  >  0 ,c2  >  O,0o  =  -oo,  and  $i  =  0, 


L  =  ci  (l  -f-  c2)  ^  f  ec,u  du 

J  —OO 

_  ci(l+c2)* 


If  the  searcher  moves  at  constant  speed  v,  then  the  right  side  of  Equation  (3-2)  must  equal  vt,  where 
t  is  the  time  required  to  go  from  the  origin  to  the  point  (p,0).  That  is,  we  must  have 


ci  (1  +  c2)* 


eC3S  =  vt. 


Now.  of  course,  p(6)  =  cieC;>9;  so  that,  considering  the  search  width  W,  the  area  of  coverage  up  to  the 
time  that  a  searcher  meets  the  boundary  of  a  circular  disk  r  <  r0  (including  recoverage  if  it  exists) 
is  just 

LW=  (1  +c|)*roW7c2.  (3-4) 

In  particular,  when  t  =  T,  one  requires  that  r0  =  R(T),  so  that  the  second  member  of  Equation  (3-4) 
must  be  equated  to  WvT.  Then  we  have 

vT=  (1  +  cD*  R(T)/c2,  (3-5) 

so  that,  necessarily,  using  the  definition  of  R{T), 


v2T2  =  (l  +  c2)  R}{T)/c\ 

=  2{l  +  cl)<r{aWvT!*)l/cl. 


Defining  c3  =  1  [c\,  one  can  readily  solve  Equation  (3-6)  for  c3  (and  hence  for  c2)  as  a  function  of 
v,T,W,  and  a.  Note,  from  Equation  (3-5),  that,  when  vT  and  R(T)  are  close  to  each  other,  c2  is 
very  large,  implying  that  the  spiral  moves  out  from  the  origin  very  rapidly  in  order  to  meet  the 
requirement  that  the  circle  p  =  R(T)  be  intercepted  at  time  T.  On  the  other  hand,  when  vT  is  large 
relative  to  R(T),  Equation  (3-5)  implies  that  c2  will  be  close  to  zero,  so  that  one  builds  up  effort  in 
the  neighborhood  of  0  before  he  transits  to  the  boundary  of  the  disk.  Observe  also  that  R(T)  <  vT  if 
and  only  if  t/3  >  4 aW<r2/nT3.  In  particular,  one  value  of  v  could  be  used  for  all  budget  times  T  >  T0, 
say. 
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It  appears  now  that  the  parameter  ci  serves  solely  as  a  scaling  parameter,  not  essentially 
affecting  our  analysis.  Therefore,  let  ct  =  1  and  c2  =  c.  We  would  like  to  examine  the  nature  of 
the  effort  function  £(r),0  <  r  <  R(T),  as  given  by  Equation  (2-27).  Note  first  that  E{ 0)  =  E'(0 )  =  0. 
Also,  E'(r )  =  nr  (i?2  -  r2)  /aa2,  which  is  positive  for  r  <  R  and  zero  for  r  =  R.  Furthermore,  one  has 
E"(r)  =  7 r  (R2  -  3r2)  /aa2.  Clearly,  then,  on  0  <  r  <  R(T),  E"(r)  has  precisely  one  root;in  addition, 
E"(r)  goes  from  positive  to  negative  values  through  the  root.  That  is  to  say,  E{r)  has  a  point  of 
inflection.  Define  £i(r;/fc)  =  kWr,k  >  0,  and  let  us  find  that  value  or  those  values  of  k  which  minimize 
the  following  cost  functional: 


Ci(*)  = 


fX'n  , 

/  [E{r)-E^r-k))  dr. 

JO 


(3-7) 


From  Equation  (3-7),  we  have 


,K(T) 

C[(k)  =  —2W  /  r[£(r)  —  kWr]  dr, 

Jo 


(3-8) 


which  vanishes  when 


k'  ~  WRXrj  J<,  rEMdr 

3  /H(T>  /  Trr3  irr5  \  J 

~  W R3(T)  J0  \2aa2  (  )  4 atr*)  dr 

ttR3(T) 

~  4 aWa2  ' 


(3-9) 


It  is  clear  from  Equation  (3-8)  (derivative  passes  from  a  negative  value  through  zero  to  a  positive 
value)  that  k *  is  a  (unique)  minimizing  point  for  the  right  side  of  Equation  (3-7).  Observe  from 
Equation  (3-9)  that  k*  >  1  whenever  -kR3{T)  >  AaWa2.  However,  by  definition  of  R(T),  we  have  that 


vT  *R3(T) 
R(T)  ~  4 aWa2 


so  that,  indeed,  Jfc*  >  1  when  R{T)  exceeds  vT ,  a  condition  which  one  naturally  requires  in  order 
for  the  searcher  to  meet  the  boundary  of  the  disk  p(t)  <  R{T)  at  time  T.  Such  being  the  case, 
k *  =  (1  +  c*)Vc*>  where  c*  is  the  solution  to  Equation  (3-5).  In  other  words,  the  best  linear  least 
squares  fit  to  the  optimal  cumulative  effort  function  E(r),  0  <  r  <  R(T),  must  be  the  logarithmic 
spiral  p(t)  =  exp(c*0(t))  which  matches  the  problem  conditions.  The  situation  is  depicted  in  Figure 
3-1. 

The  form  of  p(t)  and  0(t)  can  now  be  readily  determined  through  the  condition  p2(t)+p2(t)62(t)  = 
v2.  For  any  c  and  p{t)  =  ee,(<\  one  has 


c2e2ct02(t)  +  e2ct02(t)  =  v 2, 


so  that 


02(t) 


v2e~2ctW 
1  +  c2 


(3-10) 


3-3 


NAVSWC  TR  87-306 


FIGURE  3-1.  CUMULATIVE  EFFORT  VERSUS  OPTIMAL  LINEAR  EFFORT 
Since  9  is  to  be  a  strictly  increasing  function  of  time,  Equation  (3-10)  implies  that 

v--c8(t) 

9(t)  =  — - r.  (3-11) 

(l  +  c»)* 

Using  the  initial  condition  0(0)  =  -oo  and  separating  variables  in  Equation  (3-11),  one  finds  that 


m  =  -inf  - cvt  x- 
c  V(l+c2)’ 


(3-12) 


When  c  =  c*,  so  that  c*/  (l  +  c*2) "  =  R(T)/vT,  one  has,  from  Equation  (3-12), 

*«>  =  £l  »(Sp). 

Now  p(t)  =  c*ec*#(<)0(<),  so  that,  appealing  to  Equation  (3-13)  and  the  condition  p(T)  =  R(T),  we  have 


p(t)  = 


R(T)t 


(3-14) 


The  error  functional  Ci(ifc),  as  exhibited  in  Equation  (3-7),  is  analogous  to  that  of  Equation 
(2-26),  the  difference  being  that  integration  in  the  former  is  with  respect  to  r,  and  integration  in 
the  latter  is  with  respect  to  t.  We  can,  of  course,  make  a  change  of  variable  from  r  to  t  for  a  given 
choice  r  =  r(t).  In  particular,  suppose  that  we  use  Equation  (3-14).  Then  we  have 


W)  =  ^W(O)  -  Wvt]2  dt 


=  ^C(  r(t)), 

so  that  C(r(t))  for  the  logarithmic  spiral  chosen  is  just 

C( 


(3-15) 


(3-16) 


3-4 


NAVSWC  TR  87-306 


We  next  examine  the  behavior  of  Ci(ifc*)  as  T  tends  to  infinity.  It  will  be  shown  that  Ci(fc*)  = 
O(T^),  so  that,  from  Equation  (3-16),  C(r(i))  =  ©(T3),^)  being  given  by  the  logarithmic  spiral. 
One  has 


xfl3(r)r~|2 

4aa2 


dr 


*R3(T) 

2aa2 


*2R6(T) 
16a2  <r4 


R3(T) 

3 


(3-17) 


Let  us  consider  the  right  member  of  Equation  (3-17).  It  is 


4aV4 


r*KT)  / 

l  r'{ 


R\T)  -  R7(T)r 2  +  — 


7 r2R9{T) 
48a2er4  ' 


(3-18) 


which  reduces  to  one  term ,  namely, 


Now  R(T)  =  (2c)r(aWvT/i r)l,  so  that 


Ci(n 


x2R9(T) 
2520 a2a4' 


(3-19) 


c  ,r)  W<r)taHWvT)t 

U  '  315x4 


From  Equation  (3-16),  one  sees  that 


(3-20) 


(3-21) 


Let  us  now  compare  the  right  side  of  Equation  (3-21)  with  the  cost  obtained  for  the  optimal 
solution.  The  latter  is  the  sum  of  two  quantities  obtained  in  Chapter  2,  namely,  the  second  member 
of  Equation  (2-45)  for  J\  and  the  second  member  of  Equation  (2-50)  for  J2.  First,  let  us  consider 
J\.  Ji  depends  critically  on  the  nature  of  ti,  so  we  want  to  study  the  asymptotic  behavior  of  t\  as  T 
tends  to  infinity.  Now  is  given  by  the  right  side  of  Equation  (2-43),  and  we  observe  that  it 
depends  on  the  quantity  K  =  R(T)/vT  =  (2a)  (aW / xu3!13)^ .  Since  K  tends  to  zero  as  T  tends  to 
infinity,  suppose  that  we  let  x  =  ZKV&/8  and  study  Cos-1(-x).  Let  a  be  such  that 


Cos~\-x)=^  +  a.  (3-22) 

From  Equation  (3-22),  it  follows  that 

-x  =  cos  +  a  j  =  -  sin  a,  (3-23) 

so  that  x  =  sin  a.  Since  x  is  close  to  zero  when  T  is  large,  we  see  from  Equation  (3-22)  that  a  will 
also  be  small.  Thus,  in  Equation  (3-23),  we  may  replace  sin  a  by  a  itself  and  claim  that,  for  small  x, 

Cos~x(-x)  ~  ^  +  x,  (3-24) 


3-5 


NAVSWC  TR  87-306 


the  right  side  being  the  linear  part  of  the  Maclaurin  series  for  Cos_1(-x).  It  follows  that 

2 KTy/6  [l  (  n  ZKy/6\  An 

'■  — 3 — 008  [3  (j+-r-J  +  T 

‘iktS  ( 3*  KVe\ 

3  \2  8  ) 

2 KTV&  .  (  KsJ 6\ 

Now  sin  ~  so  that,  from  Equation  (3-25), 


(3-25) 


2 

R?(T) 

2v2T 


(  aW 


(3-26) 


Therefore,  behaves  like  T~^  as  T  becomes  large.  From  Equation  (2-45),  we  see  that  each  term  of 
Ji  tends  to  zero  as  ti  tends  to  infinity,  so  that  Jx  itself  tends  to  zero.  Indeed,  its  behavior  is  of  the 
order  of  the  last  three  terms,  so  that  Ji  goes  to  zero  like  T'l  as  T  goes  to  infinity. 

Next  consider  J2,  given  by  the  right  side  of  Equation  (2-50).  We  examine  t2,  as  provided  by 
the  second  member  of  Equation  (2-44).  Again,  t2  depends  upon  K.  Note  first  that 


Cos~l(-x)  =  n  -  Cos  '(x), 


(3-27) 


valid  for  all  x  if  we  agree  that  Cos-1  is  to  be  the  principal  arccos  operator.  Now  let  x  =  1  -  y2,  and 
consider  Cos-1(l  -  y2).  We  find  that,  for  small  positive  y  (x  close  to  l), 


Cos  *(1  —  y2)  =  Sin  1  (y\/2  -  y2) 

~  yV 2, 

where  Sin-1  is  understood  as  the  principal  arcsin  function.  Thus  we  find  that 


h  ~  T  1  — 


1  +  cos 


1  (  3\/6 k\  2n 

A'-—  ]  +  T 


^  ,  4A  ( ,  (  y/6K 

=T  ,+cos  r-  — 


-  -f  — 


8  K  .  ,  %/s K 


=r  -r 


(3-28) 


(3-29) 
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so  that  t2  tends  to  infinity  like  T  itself  as  T  tends  to  infinity.  Thus,  t3  tends  to  infinity  like  T3  &s  T 
goes  to  infinity.  Now  let  us  study  carefully  the  expansion  of  J2,  as  given  by  Equation  (2-50).  One 
finds,  after  a  simple  calculation  and  upon  making  use  of  Equation  (3-29)  and  the  definition  of  ft(T), 
that  the  first  three  terms  may  be  combined  to  give 


(T  — t2)  (WvTK7)2 
48 


(3-30) 


which  is  O  The  remaining  terms  in  the  expansion  will  likewise  be  seen  to  be  O  (t~^.  The 

result  is  that  J2,  along  with  Jt,  is  O  Therefore,  the  cost  for  the  optimal  solution  is  of  order 

T~l,  compared  to  T3  for  the  logarithmic  spiral.  It  appears  that  the  optimal  solution  is  clearly 
superior  to  the  spiral. 

Next  let  us  consider  the  elliptical  Gaussian  case.  We  want  to  define  a  curve  called  an  elliptical 
logarithmic  spiral,  which,  in  the  circular  case,  would  become  our  conventional  spiral.  As  indicated 
in  Figure  3-1,  the  radial  part  of  the  circular  spiral  was  linear;  indeed,  we  found  that  p(t)  =  R{T)t/T. 
Therefore,  for  the  elliptic  case,  again  set  p(t)  =  R(T)t/T ,  but  with  p7(t)  =  {x\{t)/<r\  +  x\{t)l<rl).  That 
is,  p[t)  is  now  the  elliptic  norm  of  r(<)  =  (xi(t),x2(t)).  At  this  point,  we  invoke  the  procedure  that  we 
employed  in  Chapter  2  to  obtain  the  optimal  search  path,  but  we  employ  that  method  in  order  to 
generate  an  elliptical  spiral.  Thus,  let  yi(t)  =  xi(t)/<T\  and  y2(t)  =  x2(t)/<r2.  In  Equation  (2-77)  replace 
f(t)  by  p2(t),  where  p(t)  =  R{T)t/T.  In  the  nonphysical  plane,  one  again  has  Equation  (2-80)  and 
Equation  (2-81),  and  we  substitute  fy(<)  =  /(<)/2ry(i)  into  Equation  (2-81),  using  our  new  f(t).  One 
then  solves  the  nonlinear  differential  equation  (2-82),  subject  to  a  convenient  initial  condition  on 
0y.  By  analogy  with  the  circular  case,  one  could  make  0y(O)  large  and  negative ,  for  example.  Also, 
one  should  choose  the  plus  sign  in  Equation  (2-82)  to  insure  that  the  curve  spirals  outward  from 
the  origin. 

Having  defined  our  elliptic  logarithmic  spiral,  let  us,  as  in  the  circular  Gaussian  case,  compare 
the  cost  function  for  the  spiral  with  the  cost  function  for  the  optimal  path.  Appealing  to  Equation 
(2-66)  and  Equation  (2-67),  we  find  that  the  cost  for  the  logarithmic  spiral  is  the  same  as  that  for 
the  circular  one,  i.e.,  we  again  have  Equation  (3-21).  Furthermore,  the  forms  for  Jx  and  J2,  together 
with  those  for  both  tx  and  t2,  are  still  valid.  Therefore,  we  find  the  same  essential  behavior  that  we 
obtained  in  the  circular  Gaussian  case. 

There  is  one  situation  in  which  the  asymptotic  behavior  of  the  cost  function  for  the  spiral 
mimics  that  for  the  optimal  solution.  Suppose  that  one  decides,  for  some  reason,  to  let  v  vary  with 
the  budget  time  T  such  that  v(T)  =  cT~a,  where  c  is  some  positive  constant.  For  the  circular  case, 
one  knows  that  R(T)  <  vT  when  v3  >  AaW<r2/irT3,  i.e.,  when  v  >  ^  /T.  This  inequality  can  be 

obtained  in  two  reasonable  ways,  namely,  (1)  by  choosing  any  constant  c  >  (4 aWo7/*)^  and  letting 
v  =  c/T  or  (2)  by  setting  v  =  (AaWcr7 /*■)$ /T° ,  0  <  a  <  1.  The  former  procedure  will  force  the  distance 
covered  in  time  T  to  be  a  constant  for  all  T  and  will  therefore  clearly  not  lead  to  a  useful  search 
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strategy.  In  contrast,  Procedure  (2)  allows  the  search  to  expand  as  the  budget  time  T  increases, 
and  we  might  like  to  adjust  a  so  as  to  minimize  the  cost.  Therefore,  let  us  assume  that 


_  ( 4aWv7\  *  ,rr* 

V"  *  ) 


0  <  c  <  1. 


(2-51) 


Then,  appealing  to  Equation  (3-21),  we  find  that 


C(r(t))  = 


2W2T3~2a 

315 


0  <  a  <  1. 


(3-32) 


Thus,  if  we  let  a  =  1  -  c,  where  c  is  small  and  positive,  we  can  almost  make  the  cost  C  linear  in  the 
time  T  and,  at  the  same  time,  assure  ourselves  that  the  search  expands  in  time.  On  the  other  hand, 
examine  the  optimal  solution  for  the  circular  case.  Consider  Ju  as  given  by  Equation  (2-45).  Noting 
that  ti  behaves  like  7<3a_1)/2  as  T  tends  to  infinity,  it  follows  that  t >8i?  =  O  (T<lla~9^2),  which,  for 
a  =  1  -e,  is  almost  0(T).  The  other  terms  in  Ji  are  all  found  to  manifest  the  same  behavior  pattern, 
so  that, loosely  speaking,  Ji  =  0{T).  Similarly,  one  finds  that  the  right  side  of  Equation  (2-50)  is 
essentially  0(T ),  implying  that  the  cost  for  the  optimal  solution  is  essentially  O(T),  provided  that 
v(T)  =  Ofl/T^l.a  =  1  -  e,  and  e  is  small  and  positive.  It  is  clear  that,  for  the  elliptic  case,  the  same 
statements  must  hold,  since,  apart  from  a  few  constants,  the  forms  of  Ju  J2.t1.t2,  and  the  cost  for 
the  spiral  +he  same.  The  choice  of  v(T )  =  0(1/T“),  with  a  close  to,  but  less  than,  1  has  the 
effect  of  eq-uilil rating  the  terms  in  Jx  and  Ji,  so  that  they  all  behave  like  T.  The  cost  per  unit  time, 
.namely,  C(r(t))/T,  is  almost  0(1). 
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CHAPTER  4 

AN  OPTIMAL  STOCHASTIC  PROCESS  FOR  THE  CASE 
OF  A  MOVING  TARGET 

The  main  effort  of  this  report  has  been  devoted  to  ascertaining  the  nature  of  optimal  paths 
for  acquiring  targets  known  to  be  stationary  in  some  two-dimensional  Euclidean  reference  frame. 
If  the  target  is  moving,  it  is  possible  to  adapt  the  theory  we  have  developed  so  far  in  order  to 
maximize  one’s  chance  of  finding  it  (or,  equivalently,  to  minimize  the  mean  time  to  acquire  it).  In 
fact,  we  shall  be  led  to  a  stochastic  process,  as  opposed  to  a  deterministic  one,  which  one  or  more 
search  vehicles  must  follow  within  the  time  T  allowed. 

The  idea  is  conceptually  very  simple.  Suppose,  as  shown  in  Figure  4-1,  that  we  have  located 
the  nominal  target  position  at  some  point  in  time,  which  we  shall  denote  by  t=0.  That  is  to  say, 
we  know  the  nominal  targ-'t  range  RT  and  bearing  Bt  at  t  =  0.  If  we  knew  the  target’s  velocity 
vector  vT  and  were  willing  to  travel  at  a  closing  speed  v0  to  intercept  the  most  likely  (nominal) 
position  of  the  target  after  </  =  RT/v0  units  of  time,  then  we  would  have  a  deterministic  trajectory 
for  acquisition.  In  fact,  we  see  that 


xm(<)  =  (vT  +  v0)</  +  vT(<  -  </)  +  x,(t  -  ti) 

=  yr0tj +  \Tt  +  x,(t-ij),  0  <t-tj<T,  (4-1) 

where  xm(t)  is  our  path  for  acquiring  the  moving  target,  tj  is  nominal  point  intercept  time,  v0  is  own 
ship  (or  possibly  own  missile)  velocity  when  vT  =  0,  vT  is  target  velocity,  and  x,(<-</)  is  our  trajectory 
plan  for  detecting  a  stationary  target.  Since  our  search  in  the  moving  frame  of  reference  (the  frame 
stationary  relative  to  the  target)  must  start  at  t  =  tj ,  the  total  time  budget  T  is  calculated  relative 
to  such  intercept  time.  Shown  in  Figure  4-1  are  some  of  the  level  curves  of  the  target’s  Gaussian 
location  distribution  in  the  moving  reference  frame.  The  searcher’s  total  velocity  is  v0  +  vt- 

So  far  there  is  only  one  distribution  function  for  the  target  which  needs  to  be  considered, 
namely,  that  of  the  target  within  its  own  frame  of  reference.  Since  we  know  the  target’s  velocity,  we 
need  only  travel  with  the  target  while  we  move  toward  the  intercept  point  (mode  of  the  Gaussian 
distribution)  and  continue  to  follow  the  target  motion  thenceforward. 

Now  let  us  suppose  that  vt  is  not  known  with  certainty,  but  that  it  is  some  constant  vector, 
whose  distribution  function  is  known  a  priori.  One  example  of  such  a  distribution  function  would  be 
the  following:  The  target  speed  vt  lies  between  two  known  limits  vl  and  vy,  and  the  target  heading 
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TARGET  DISTRIBUTION 
LEVEL  CURVE 

NOMINAL  TARGET 
_ POINT  (t=  0) 


FIGURE  4-1.  SEARCHER-TARGET  GEOMETRY  FOR  MOVING 
TARGET  DETECTION  PROBLEM 

Or  is  between  9L  and  9u ■  Furthermore, 

p(vt,9t)  =  7 - T72 - T— r,  vL  <  vt  <  vv  ,9 L  <  9T  <  9V .  (4-2) 

{vu  —  vl){Qu  ~0l) 

That  is,  the  target  velocity  is  uniform  over  a  known  range  of  angles  and  speeds.  Then  we  may 
conceive  of  Equation  (4-1)  as  a  stochastic  path  process  for  acquiring  the  target,  and  we  could,  in 
fact,  design  strategies  for  target  acquisition  based  on  our  knowledge  of  the  process.  If,  for  example, 
one  is  limited  to  a  single  search  vehicle,  then  it  might  be  reasonable  to  follow  the  mean  of  the 
process  and  hope  for  the  best.  Thus, 

xm(t)  =  £(xm(<))  =  x0l,  +  tE(vT)  +  x,(<  -  </).  ti  <t  <tj  +  T,  (4-3) 

so  that  xm(!)  is  again  a  deterministic  search  mechanism.  On  the  other  hand,  if  several  vehicles 
could  be  released  from  point  O  in  Figure  4-1  at  roughly  the  same  time,  we  could  try  to  follow  the 
stochastic  process  itself  at  some  level.  That  is,  let 

xmHt)  =  vo</  +  h  +  x,(<  -  </),  1  <  k  <  NV,  (4-4) 

where  Xm\t)  is  realization  k  of  stochastic  motion  xm(i)  for  search  vehicle  E,  NV  being  the  number  of 
vehicles  available. 
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One  might  imagine  several  variants  of  the  stochastic  motion  we  have  presented  here,  depending 
upon  known  characteristics  of  target  motion  and/or  searcher  motion.  For  example,  suppose  that 
the  target  is  known  to  be  decelerating  over  time  in  some  given  direction,  having  started  from  either 
a  known  or  unknown  speed  at  time  0.  Then  the  generic  equation  (4-1)  is  still  operative;  only  the 
distribution  given  by  Equation  (4-2)  need  be  altered.  One  candidate  might  be 


p(v^\eT,dT)  = 


(4-5) 


\vu  ~  v‘L){0u  ~  Oi.){du  ~  d/.)  ’ 

where  v is  the  initial  target  speed  and  dT  is  the  constant  target  deceleration  (a  negative  deceleration 
is  interpretable  as  an  acceleration). 

Of  course,  one  may  learn  over  time  (if  there  is  enough  time  to  do  so),  thus  allowing  one 
to  distribute  a  number  of  search  vehicles  on  a  tentative  basis.  These  vehicles  might  be  cheap, 
dispensable  items.  Once  the  target’s  motion  is  known  with  greater  certainty,  a  more  expensive 
search  mechanism  could  be  employed,  perb?ns  pfong  the  expected  course  (presumably  then  having 
small  variance  characteristics). 
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CHAPTER  5 

SUMMARY  AND  FUTURE  WORK 

We  have  shown  in  this  report  how  one  can  use  optimal  control  theory  to  design  tracking 
procedures  which  maximize  the  acquisition  probability  for  either  stationary  or  moving  targets.  It 
is  assumed  that  one  is  given  a  time  budget  T  in  which  to  operate,  that  one  may  choose  his  search 
speed  v  within  reasonable  operating  limits,  and  that  one  will  initiate  search  at  the  most  likely  target 
point.  We  have  compared  the  efficiency  of  optimal  search  plans  with  logarithmic  spiral  plans,  which 
are  obviously  suboptimal,  and  we  found  the  following:  If  speed  is  not  chosen  carefully  as  a  function 
of  budget  time,  the  optimal  search  is  clearly  superior  to  the  spiral  search.  However,  if  one  is  careful 
to  adjust  the  search  speed  and,  in  particular,  if  one  moves  relatively  slowly,  both  optimal  search 
and  logarithmic  spiral  search  have  similar  characteristics  for  large  budget  times. 

If  one  does  not(or  cannot)  initiate  search  from  the  most  likely  target  point,  he  could  be  forced 
to  pursue  a  different  track.  However,  in  certain  circumstances  it  would  be  expected  that  a  simple 
reversal  of  an  optimal  path  mechanism  or  a  logarithmic  spiral  search  would  be  applicable.  For 
example,  in  a  radially  symmetric  situation,  as  would  be  encountered  when  the  targt.  distribution 
is  circular  normal,  if  one  had  to  start  search  on  the  boundary  of  some  disk  at  a  prescribed  point , 
working  into  the  disk,  he  should  probably  reverse  the  path  which  would  have  normally  evolved  from 
the  center  (where  the  target  was  most  likely  to  be).  On  the  other  hand,  an  elliptical  distribution 
would  present  a  different  situation  if  one  had  to  start  at  a  point  not  on  the  minor  axis.  That  case 
has  not  been  studied  in  this  report. 

Future  work  would  address  the  important  issue  of  tracking  targets  in  the  presence  of  either 
false  targets  or  aggressor  targets  whose  influence  regions  overlap  those  of  the  target  of  interest. 
Also,  numerical  comparisons  between  optimal  and  suboptimal  paths  developed  here  might  be  of 
interest,  as  well  as  comparisons  between  them  and  paths  actually  utilized  in  practical  situations. 
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